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1.  INTRODUCTION 


This  user  guide  explains  the  new  CHIEF2000  computer  program.  The  program  was 
tested  on  machines  running  Windows  95/98  ®,  Windows  NT  ® ,  and  DEC  OSF1.  The  CHIEF 
program  was  originally  developed  in  the  1960s  to  compute  the  acoustic  radiation  from  an 
arbitrary  shaped  radiating  body.  The  program  was  extensively  updated  in  1988  and  a  users 
manual  was  published  [Benthien,  Barach,  and  Gillette,  1988].  This  CHIEF  manual  is  still 
an  excellent  reference  for  those  interested  in  the  details  of  the  CHIEF  program  and  the 
underlying  theory.  The  CHIEF  manual  was  included  in  the  recently  released  CD-ROM 
document  [Benthien  and  Barach,  1998].  The  appendix  summarizes  CHIEF  program  theory. 

The  original  CHIEF  program  and  the  1988  update  were  essentially  a  library  of  Fortran 
subroutines  that  were  called  from  a  user-generated  control  routine.  This  approach  allowed 
the  user  more  flexibility,  but  required  a  Fortran  compiler  to  compile  and  link  the  control 
routine  with  the  library  routines.  The  presence  of  a  Fortran  compiler  was  not  a  big  problem 
on  most  main  frame  computers.  However,  with  the  increasing  use  of  PCs  and  workstations 
for  numerical  computation,  the  presence  of  a  Fortran  compiler  could  no  longer  be  assumed. 
Thus,  the  new  CHIEF  2000  program  is  a  compiled  executable  program  that  relies  on  a  user¬ 
generated  ASCII  text  file  for  input.  While  this  approach  gives  up  some  flexibility,  we  have 
included  as  many  features  of  the  previous  versions  of  CHIEF  as  possible. 

This  user  guide  describes  the  setup  of  the  user-generated  input  file.  Section  2  describes 
the  use  of  the  program  and  the  general  organization  of  the  input  file.  Sections  3-5  describe 
the  main  subdivisions  of  the  input  file.  Section  3  describes  the  parameter  section  in  which 
the  user  can  define  parameters  for  use  in  the  remaining  sections  of  the  input  file.  Section 
4  describes  the  geometry  input  section  in  which  the  user  specifies  the  radiating  surface 
and  its  subdivision.  Section  5  describes  the  control  section  in  which  the  user  specifies  the 
quantities  to  be  computed  and  the  form  of  the  output.  Section  6  describes  the  format  of 
the  various  possible  output  files.  Section  7  contains  many  example  problems  along  with  the 
corresponding  input  files. 

Much  of  the  material  in  this  user-guide  is  reference  material;  therefore  this  document  is 
not  designed  to  be  read  from  start  to  finish  as  one  would  read  a  novel.  The  authors  feel 
that  the  program  can  be  best  learned  by  studying  the  example  problems  near  the  end  of  the 
manual.  A  suggested  reading  plan  is  given  below: 

1.  Read  Introduction. 

2.  Read  General  Organization  section. 

3.  Read  Parameter  section. 

4.  Read  introductory  portions  of  the  Geometry  and  Control  sections. 

5.  Work  through  example  problems. 

6.  Refer  to  the  command  descriptions  in  the  Geometry  and  the  Control  sections  as  needed. 
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2.  GENERAL  ORGANIZATION 


The  program  can  be  run  by  typing  the  command 
chief  <input  file  name>.inp 

at  the  command  prompt.  In  Windows  95/98®  or  Windows  NT®,  the  user  should  go  to  an 
MS-DOS  ®  window  to  enter  the  command.  The  input  file  name  can  be  any  legal  DOS  file 
name  (maximum  of  eight  characters  with  no  spaces).  The  .inp  extension  can  be  omitted  in 
the  command  line.  If  the  run  completes  normally,  the  message  Normal  Termination  will 
appear  on  the  screen. 

The  input  file  (cinput  file  name>.inp)  is  an  ASCII  text  file  generated  by  the  user.  It 
must  not  contain  any  control  characters.  The  file  is  divided  into  the  following  three  major 
sections: 

1.  Parameter  definition 

2.  Geometry  specification 

3.  Control  logic 

In  the  parameter  section,  the  user  can  assign  parameter  names  to  values  required  for 
input  in  the  following  sections.  The  parameter  definitions  can  involve  arithmetic  expressions 
as  well  as  previously  defined  parameters.  It  is  not  necessary  to  define  any  parameters  since  all 
required  input  values  can  be  either  parameters  or  numerical  values.  The  parameter  section  is 
provided  for  the  user’s  convenience.  In  the  geometry  section,  the  user  specifies  the  radiating 
surface,  the  subdivision  of  the  surface,  and  any  symmetry  conditions.  In  the  control  logic 
section,  the  user  can  specify  the  quantities  to  be  computed,  provide  necessary  input  data, 
and  specify  the  type  of  output.  A  sample  input  file  is  shown  below. 
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!  Parameter  section 


nblks=8 

rotlim=pi/nblks 

radius=l .2e-l  !  radius  of  sphere 

freql=100 

frequ=200 

freqinc=10 

end*parameters 

!  Geometry  section 

rotational , nblks 

sphere , radius 

0 , pi , -r ot 1 im , rot 1 im 

12,1,4,16 

3, 0,0,0 

end*geometry 

!  Control  logic  section 

frequency , f reql , f requ , f reqinc 

fluid, 995, 1489 

rhs  symmetry 

velocities , sample . vel 

surface  pressures , 1 

far-f ield, 1,3 

0,90,5,0,0,0 


A  detailed  explanation  of  each  input-file  section  will  be  given  later.  Note  that  the  pa.ra.m- 
eter  definition  section  must  be  terminated  with  a  statement  beginning  with  end*p  and  the 
geometry  section  must  be  terminated  with  a  statement  beginning  with  end*g.  The  control 
logic  section  is  terminated  by  an  end  of  file.  Blank  lines  are  ignored  in  the  input  file,  as  are 
any  comments.  Comments  must  begin  with  an  exclamation  point  (!).  All  commands  and 
parameters  are  case  insensitive.  Thus,  NBLKS  is  the  same  parameter  as  nblks,  and  Sphere 
is  the  same  command  as  sphere.  It  is  also  important  that  consistent  units  be  used  for  all 
input  quantities.  The  program  by  default  takes  the  fluid  density  and  sound  speed  to  be  1000 
Kg/m3  and  1500  m/  sec.  If  the  user  desires  to  use  a  system  of  units  other  than  mks,  then 
these  quatities  must  be  changed  to  values  appropriate  for  that  system  of  units.  Data  read 
in  from  files  must  also  be  in  the  chosen  consistent  set  of  units. 
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Unlike  the  original  version  of  CHIEF,  this  new  version  has  very  few  predetermined  size 
limits.  Most  arrays  used  by  the  program  are  allocated  dynamically  within  the  program 
based  on  inputs  in  the  input  file.  There  are  a  few  arrays,  however,  that  are  still  statically 
dimensioned.  Thus,  the  maximum  size  of  these  arrays  cannot  be  exceeded  without  modifying 
the  program  and  recompiling.  Table  1  shows  these  limits.  We  do  not  feel  that  these  limits 
should  impose  any  serious  limitations. 


Table  1.  Array  size  limits  in  CHIEF2000. 


Maximum  finite  element  node  number 

5000 

Maximum  number  of  interior  points 

100 

Maximum  number  of  point  sources  in  point  source  check 

100 

Maximum  number  of  Gaussian  quadrature  points  on  a  sur¬ 
face  subdivision  in  either  of  the  two  surface  coordinate  di¬ 
rections 

64 

Maximum  total  number  of  Gaussian  quadrature  points  on  a 
surface  subdivision 

512 

Maximum  order  of  rotational  symmetry 

120 

In  addition  to  any  user-specified  output  files,  the  program  always  creates  a  file  with  the 
name  <input  file  name> .  out  that  contains  an  echo  of  all  input  file  lines,  error  diagnostics, 
and  user-specified  output.  In  the  parameter  and  geometry  sections,  the  program  finds  as 
many  input  errors  as  it  can  before  stopping.  However,  the  program  does  not  go  on  to  the 
next  section  if  there  are  any  errors  found  in  the  current  section.  In  the  control  section,  the 
program  terminates  if  any  input  error  is  found.  Termination  helps  prevent  calculations  from 
being  performed  with  inconsistent  or  missing  input. 


The  program  can  generate  many  output  files.  These  files  all  have  the  same  root  file 
name  as  the  input  file.  Table  2  lists  the  extensions  of  these  output  files  along  with  a  brief 
description  of  their  contents. 


Table  2.  Possible  output  files. 


File  Extension 

File  Type 

File  Contents 

.sps 

ASCII 

Surface  pressures  on  each  subdivision 

.ffp 

ASCII 

Far-field  angles  and  pressures 

.nfp 

ASCII 

Near-field  coordinates  and  pressures 

.zmx 

Binary 

Radiation  impedance  matrix 

.pmf 

Binary 

Fax-field  pattern  matrix 

.pmn 

Binary 

Near-field  pattern  matrix 

.abs 

Binary 

A  and  B  surface  matrices 

.abf 

Binary 

A  and  B  far-field  matrices 

.abn 

Binary 

A  and  B  near-field  matrices 
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3.  PARAMETER  SECTION 


In  the  parameter  section,  the  user  can  define  parameters  that  can  be  used  elsewhere  in 
the  CHIEF  ASCII  input  file  in  place  of  numerical  inputs.  Only  a  single  parameter  may  be 
defined  on  each  line  of  the  parameter  section,  and  the  parameter  definition  must  have  the 
form 

<parameter  name>=<expression> ♦ 

The  parameter  name  may  be  any  sequence  of  alpha-numeric  characters  or  underscore 
characters  (_  )  subject  to  the  following  restrictions: 

•  The  first  character  must  be  a  letter. 

•  The  length  of  the  parameter  name  can  not  exceed  24  characters. 

•  The  parameter  name  cannot  contain  any  spaces. 

Case  is  not  important  in  parameter  names  (e.g.,  Nblks  is  the  same  as  nblks). 

All  parameter  values  are  interpreted  by  the  program  as  real  numbers.  However,  the 
nearest  integer  is  used  in  places  where  an  integer  input  is  expected.  The  user  may  insert 
comments  on  a  separate  line  or  at  the  end  of  a  line  by  preceding  the  comment  with  an 
exclamation  point  (!). 

The  expression  on  the  right-hand  side  is  subject  to  the  following  conditions: 

•  The  four  arithmetic  operations  (+,  — ,  *,  /)  and  exponentiation  (A)  are  allowed. 

•  The  sine  or  cosine  of  an  arithmetic  expression  is  allowed  (i.e. ,  sin(x"2)  or  cos(x+y)). 

•  The  operands  can  be  numerical  values  or  previously  defined  parameters. 

•  Parenthetical  grouping  (...)  and  nested  operations  may  be  used. 

•  An  exponent  may  be  any  real  number. 

•  Scientific  notation  is  allowed  (i.e.,  1.6e-6). 

•  The  entire  fine  must  not  exceed  132  characters. 

Expressions  satisfying  the  above  conditions  are  also  allowed  in  the  geometry  and  the  control 
sections  wherever  a  numeric  input  is  expected. 

Spaces  may  be  freely  used  in  parameter  definitions  except  within  parameter  names.  The 
program  ignores  them.  There  is  one  parameter,  pi,  that  has  been  predefined  by  the  program 
and  may  be  used  anywhere  in  the  input  file,  including  the  definition  of  other  parameters. 
The  usual  rules  for  precedence  of  operations  apply,  for  example, 
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•  Operators  of  equal  precedence  are  evaluated  in  the  usual  left  to  right  order. 

•  Exponentiation  has  precedence  over  addition,  subtraction,  multiplication,  and  division. 

•  Multiplication  and  division  have  precedence  over  addition  and  subtraction. 

•  Unary  plus  and  minus  operations  (i.e.,  -1.0  or  +5.6)  have  the  same  precedence  as 
exponentiation. 

The  parameter  section  must  end  with  a  line  whose  first  five  nonblank  characters  are  end*p. 
The  following  is  an  example  of  a  valid  parameter  section. 


rho=1000 

c=1500. 


diam 

height 

lengthl 

length2 


1.5e2  !  This  is  a  comment 

4.67 

( (height* . 5) "2+ (diam/2) “2) “ . 5 
lengthl “(1/2) 


theta_upper=pi/4 
costu=cos (theta_upper) 
theta_lower=-theta_upper 


freq_upper=3e3 
f  req_lower=f  r eq_upper * 1 e- 1 


near_f ield_r=lengthl*c/freq_lower 
near _f i eld_thet  a=pi *5E- 1 


end*parameters 


The  following  example  contains  some  invalid  parameter  definitions. 

rho  f luid=1000  (space  in  name) 

c(l)=1500  (parenthesis  in  name) 

height=4.67  width=2.0  (more  than  one  parameter  on  line) 

lengthl=( (height* . 5) “2+ (diam/2) “2) “ . 5 

(parameter  used  on  the  right  side  before  its  definition) 
ldiam=1.5,  (name  begins  with  a  numeral  and  comma  illegal) 
ixmer*radius= . 5  (*  illegal  symbol  in  name) 

outer-radius=2 . 5  (-  illegal  symbol  in  name) 

end*param 

A  line  specifying  the  end  of  the  parameter  section  (i.e.,  end*params)  must  be  included 
even  if  no  parameters  are  defined. 
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4.  GEOMETRY  SECTION 


CHIEF  calculates  radiation  from  one  or  more  closed  surfaces  in  3-space,  or  from  one  or 
more  planar  regions  in  an  infinite  rigid  baffle.  In  this  section  of  the  manual,  we  describe  how 
the  user  specifies  the  geometry  of  the  radiating  surface.  If  the  surface  has  one,  two,  or  three 
planes  of  reflective  symmetry  or  a  finite  order  of  rotational  symmetry,  then  the  user  only 
needs  to  input  the  portion  of  the  surface  corresponding  to  one  symmetry  block. 

The  commands  that  specify  the  surface  geometry  may  be  given  in  any  order,  except  the 
command  specifying  the  symmetry  type  (NONE,  REFLECTIVE,  ROTATIONAL)  must  be 
the  first  command  in  the  geometry  section.  If  there  are  no  underscore  characters  in  the 
command  word,  the  program  uses  only  the  first  three  characters.  If  there  is  an  underscore 
character,  the  program  uses  the  first  letter,  the  underscore  character,  and  the  first  character 
after  the  underscore  to  specify  the  command.  Other  forms  of  the  command  words  are 
perfectly  valid  as  long  as  the  three  characters  specifying  the  command  are  the  same. 

Each  command  line  consists  of  a  command  or  a  command  followed  by  one  or  more 
arguments.  Some  command  lines  must  be  followed  by  additional  lines  of  data.  Commands, 
arguments,  and  data  items  are  separated  by  commas.  Parameters  defined  in  the  parameter 
section  and  expressions  involving  these  parameters  may  be  used  for  any  numerical  input  in 
this  section.  As  always,  blank  lines,  comments,  and  spaces  are  ignored. 

CHIEF  has  many  predefined  surfaces  such  as  spherical,  cylindrical,  ellipsoidal,  etc.  The 
radiating  surface  can  be  defined  in  terms  of  portions  of  several  of  these  predefined  surfaces. 
CHIEF  also  allows  the  surface  to  be  defined  in  terms  of  nodes  and  elements  as  in  a  finite 
element  problem.  The  user  must  supply  inputs  specifying  the  various  nonoverlapping  surface 
regions  that  make  up  the  total  radiating  surface.  In  the  case  of  finite  element  type  input, 
each  element  is  considered  as  a  separate  surface  region.  Each  surface  region  is  defined  by  a 
one-to-one  smooth  mapping  from  a  rectangular  planar  region  into  3-dimensional  space.  The 
points  on  the  surface  as  well  as  any  field  points  are  defined  in  terms  of  a  global  rectangular 
Cartesian  coordinate  system,  xx,  x2,  x3.  Points  in  the  base  rectangular  region  are  defined  in 
terms  a  pair  of  coordinates,  u,  v.  The  values  u  and  v  are  sometimes  called  surface  coordinates. 
Each  surface  region  is  defined  by  three  relations  of  the  form: 


Xx  =  Xx (u,  V )  X2  =  X2(u,  V )  Xz  —  Xz{u,  V ). 


Each  predefined  geometry  is  actually  defined  in  terms  of  a  local  x ,  y,  z  Cartesian  coordi¬ 
nate  system  that  is  related  to  the  global  coordinate  system  by  a  translation  and/or  a  cyclic 
permutation  of  the  coordinate  axes.  The  translation  is  specified  by  a  vector,  t,  from  the 
global  origin  to  the  local  origin.  The  cyclic  permutation  of  the  coordinate  axes  is  specified 
by  giving  the  axis  ( Iz )  in  the  global  system  that  corresponds  to  the  2-axis  in  the  local  system. 
The  user  input  for  this  specification  is  the  integer  Iz.  Figures  1-3  show  the  relation  between 
the  global  and  local  coordinate  systems  for  Iz  =  1, 2,  3.  Note  that  Iz  =  3  corresponds  to  no 
rotation  at  all. 
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Xi 


Figure  1.  Global  and  local  coordinate 
systems  for  Iz  —  1. 


x, 


Figure  2.  Global  and  local  coordinate 
systems  for  Iz  —  2. 


z 


Figure  3.  Global  and  local  coordinate 
systems  for  Iz  =  3. 


Thus,  each  surface  region  is  defined  by  relations  of  the  form: 

x  =  x(u,v)  y  —  y{u,v)  z  =  z(u,v), 
and  the  specification  of  Iz  and  t. 

In  addition  to  specifying  a  relation  between  the  local  and  global  coordinate  systems,  the 
input  quantity  Iz  is  also  used  to  help  specify  the  outward  normal  vector  to  the  surface  region. 
Prom  the  defining  equations  for  a  surface  region,  CHIEF  can  compute  a  normal  vector  at  each 
point  of  the  region,  but  it  cannot  determine  whether  this  normal  vector  is  inward  or  outward. 
The  user  specifies  this  by  the  sign  of  Iz.  A  plus  sign  accepts  the  choice  made  by  the  CHIEF 
program  and  a  minus  sign  reverses  the  direction  of  the  surface  normal  on  this  surface  region. 
For  example,  in  the  case  of  a  circular  planar  surface  region,  the  program  computes  the  normal 
pointing  in  the  positive  ^-direction.  This  direction  may  or  may  not  be  the  correct  choice  for 
the  surface  region  under  consideration.  If  the  radiating  surface  was  like  a  tin  can  with  axis 
of  rotation  along  the  z-axis,  then  the  outward  normal  would  be  in  the  positive  z-direction 
on  the  top,  but  in  the  negative  z-direction  on  the  bottom.  Thus,  the  user  would  specify 
a  negative  Iz  for  the  bottom  surface.  The  user  can  also  specify  a  breakup  of  each  surface 
region  into  subareas  (on  which  the  pressure  and  normal  velocity  will  be  assumed  constant) 
by  specifying  a  uniform  breakup  in  the  uv  rectangular  region.  Conventional  wisdom  says 
that  all  subdivisions  should  be  less  than  1/4  wavelength  with  something  on  the  order  of  1/8 
wavelength  desireable.  As  far  as  integration  over  the  subdivisions,  four  Gaussian  quadrature 
points  in  1/8  wavelength  seems  to  give  sufficient  accuracy. 

The  first  statement  in  the  geometry  section  must  be  one  of  the  CHIEF  symmetry  op¬ 
tions  (ROTATIONAL,  REFLECTIVE,  or  NONE)  along  with  any  associated  parameters  such  as  the 
number  of  planes  of  reflective  symmetry.  Following  the  symmetry  specification,  the  surface 
regions  may  be  defined  in  any  order.  However,  the  ordering  will  influence  the  numbering  of 
the  subdivisions.  The  numbering  begins  with  the  first  region  specified  and  continues  consec¬ 
utively  with  the  other  regions  in  the  order  specified.  Within  a  region,  the  numbering  starts 
with  the  subarea  having  the  smallest  u,  v  values.  The  numbering  continues  through  the  re¬ 
gion  in  the  direction  of  increasing  u  and  v,  with  the  variation  in  v  occuring  before  variation 
in  u.  Figure  4  shows  an  example  of  the  numbering  for  three  adjacent  planar  surface  regions. 

The  numbering  proceeds  sequentially  through  the  various  symmetry  blocks  and  the  num¬ 
bering  within  each  symmetry  follows  the  same  pattern  as  in  the  first  block.  The  numbering 
of  symmetry  blocks  for  reflective  symmetry  is  shown  in  tables  3-5. 

The  geometry  section  ends  with  a  statement  like  end*geometry,  where  only  the  first 
part,  end*g,  is  significant  to  the  program.  Commands  specifying  the  surface  geometry  are 
described  on  the  following  pages.  The  syntax  for  each  command  is  given  along  with  a 
description  of  the  parameters  and  data  items  that  accompany  the  command.  For  commands 
that  do  not  contain  an  underline  character,  the  program  only  looks  at  the  first  three  letters. 
For  commands  that  contain  an  underline  character,  the  program  looks  at  the  first  character, 
the  underline  character,  and  the  first  character  after  the  underline  character.  Alternate 
forms  of  the  commands  are  allowed  as  long  as  the  three  significant  characters  are  in  the 
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Table  5.  Three  planes  of  symmetry 
(xi  =  0,  X2  =  0,  and  x$  =  0). 


Block 

Octants 

Xi  x2  x3 

1 

+  +  + 

2 

-  +  + 

3 

-  -  + 

4 

+  -  + 

5 

+  -  - 

6 

- 

7 

-  +  - 

8 

+  +  - 

right  place.  Optional  parameters  are  enclosed  in  brackets.  For  those  commands  that  define 
a  surface  region,  the  mapping  from  the  u,  v  rectangular  region  to  the  surface  is  described. 
The  user  needs  to  understand  what  u  and  v  are  in  order  to  understand  the  numbering  of  the 
surface  subdivisions. 
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reflective  symmetry 

specifies  planes  of  reflective  symmetry 


Command: 

N 


reflective  symmetry,  N 

Number  of  planes  of  reflective  symmetry  (1,2,  or  3) 
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rotational  symmetry 

specifies  finite  order  rotational  symmetry 


Command: 

N 


rotational  symmetry,  N 

Specifies  TV-fold  rotational  symmetry 
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none 

specifies  no  symmetry 

Command:  none 


rectangular-planar 

rectangular  region  on  surface  parallel  to  the  x,  y  plane 

Command:  rectangular_planar,  z0 

Ul,  Uu,  Vl,  Vu 
Nu,  NVy  Qui  Qv 
h )  [h,  h,  ^3] 

Constant  value  of  z  in  local  coordinate  system 
Lower  and  upper  limits  for  the  surface  coordinate  u 
Lower  and  upper  limits  for  the  surface  coordinate  v 

Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 
region 

Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 
normal 

Components  of  the  translation  vector  from  the  global  origin  to  the  local 
origin  (default  is  ti  =  t2  =  <3  =  0) 

Surface  definition: 

The  local  coordinate  system  must  be  set  up  so  that  the  z-axis  is  normal  to  the  plane 
of  the  surface.  In  this  case  u  is  equal  to  x  and  v  is  equal  to  y.  The  local  coordinates  and 
normal  vector  are  related  to  the  u,  v  coordinates  by 

x  =  u  nx  —  0 

ny  =  0 

nz  =  1 


y  =  v 
z  =  z0 


circular_planar 

circular  region  on  surface  parallel  to  the  x,  y  plane 

Command:  circular_planar,  z0 

Ul,  Uu,  VL,  Vu 
Nm  Qui  Qv 

h,  [tl,  ^2)  £3] 

Constant  value  of  2  in  local  coordinate  system 
Lower  and  upper  limits  for  the  surface  coordinate  u 
Lower  and  upper  limits  for  the  surface  coordinate  v 

Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 
region 

Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 
normal 

Components  of  the  translation  vector  from  the  global  origin  to  the  local 
origin  (default  is  ti  =  t2  —  t3  =  0) 

Surface  definition: 

The  local  coordinate  system  must  be  set  up  so  that  the  2-axis  is  normal  to  the  plane 
of  the  surface.  In  this  case  u  is  the  radial  distance  from  the  2-axis  and  v  is  the  angle  from 
the  x-axis  in  the  xy  plane.  The  local  coordinates  and  normal  vector  are  related  to  the  u,  v 
coordinates  by 

x  =  u  cos  v  nx  =  0 

y  =  usinv  ny  =  0 

z  =  z0  nz~  1 
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ellipticaLplanar 

elliptical  region  on  surface  parallel  to  the  x,  y  plane 


Command:  ellipticaLplanar,  a,  b,  zq 

Ul,  Ujj,  Vl,  Vu 
Nu,  Nvi  Qu,  Qv 

Iz ,  [ti,  t2,  ta] 


a 

b 

zo 

Ul,  uu 
vl,  Vu 
■Nw  Nv 


Qu,  Q 


v 


Iz 


tl,  t2,  <3 


Half  the  length  of  the  major  axis  of  the  ellipse 
Half  the  length  of  the  minor  axis  of  the  ellipse 
Constant  value  of  z  in  local  coordinate  system 
Lower  and  upper  limits  for  the  surface  coordinate  u 
Lower  and  upper  limits  for  the  surface  coordinate  v 

Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 
region 

Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 
normal 

Components  of  the  translation  vector  from  the  global  origin  to  the  local 
origin  (default  is  ti  =  t2  =  £3  =  0) 


Surface  definition: 

The  local  coordinate  system  must  be  set  up  so  that  the  2-axis  is  normal  to  the  plane  of 
the  surface.  In  this  case  u  =  r/ /  where  r  is  the  average  of  the  radial  distances  of  a  point 
from  the  two  foci  and  /  is  half  the  interfocal  distance.  The  range  of  u  is  1  <  u  <  a/ \/a2  —  b2. 
The  family  of  curves  for  which  u  is  a  constant  consists  of  ellipses  having  the  same  foci.  The 
ellipse  corresponding  to  u  =  1  is  the  degenerate  ellipse  consisting  of  the  interval  [-f,f]  on 
the  rr-axis.  The  ellipse  corresponding  to  u  =  a/y/a 2  -  b2  is  the  boundary  ellipse  with  major 
diameter  of  2 a  and  minor  diameter  of  2b.  Thus,  the  upper  limit  uv  of  the  variable  u  is  almost 
always  a/y/a 2  -  b2.  On  an  ellipse  of  constant  u  the  variable  v  is  related  to  the  angle  6  from 
the  x-axis  in  the  xy  plane  by 

v  =  tan-1  (  .  U  tan# 

\y/tf=T 
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The  local  coordinates  and  normal  vector  are  related  to  the  u,  v  coordinates  by 

x  =  fu  cos  v  nx  =  0, 

y  =  f\/u2  -  lsinv  7^  =  0, 

z  =  nz  =  1, 

where  /  =  yfa?  —  b2. 

Equal  increments  in  v  do  not  produce  equal  increments  in  9  and  vise  versa.  However, 
quadrants  in  the  rry-plane  do  come  out  the  same,  as  can  be  seen  from  the  following  table: 


V 

X 

y 

9 

0 

fu 

0 

0 

7r/2 

0 

f\/u2  - 1 

7t/2 

7T 

-fu 

0 

7T 

3-7t/2 

0 

-fVu2  -  1 

37r/2 
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circular-cylinder 

cylindrical  surface  with  axis  along  z-axis 

Command:  circular -cylinder,  a 

Ul,  Uu,  VL,  Vxj 

Nu,  NV)  Qut  Qv 

hi  [^1>  h,  *3] 

radius  of  cylinder 

Lower  and  upper  limits  for  the  surface  coordinate  u 
Lower  and  upper  limits  for  the  surface  coordinate  v 

Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 
region 

Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 
normal 

Components  of  the  translation  vector  from  the  global  origin  to  the  local 
origin  (default  is  t\  =  £2  =  ^3  =  0) 

Surface  definition: 

The  local  coordinate  system  must  be  set  up  so  that  the  z-axis  coincides  with  the  axis  of 
the  cylinder.  In  this  case  u  —  z  and  v  is  the  angle  from  the  a>axis  in  the  xy  plane.  The  local 
coordinates  and  normal  vector  are  related  to  the  u,  y  coordinates  by 

x  =  a  cos  v 
y  =  a  sin  v 
z  =  u 


nx  =  cos  v 
ny  =  sin  v 
nz—  0 
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elliptic-cylinder 

elliptic  cylindrical  surface  with  axis  along  2-axis 


Command:  elliptic_cylinder,  a,  b 

Ul,  Uu,  Vl,  Vu 

Nu,  Nv,  Q 

Ul  Qv 

Izi  [^1>  til  ^3] 


a 

b 

Ul,  Uu 

ul ,  vu 
Nu,  Nv 


Qui  Qv 


^1)  ^2j  tz 


Half  the  major  diameter 
Half  the  minor  diameter 

Lower  and  upper  limits  for  the  surface  coordinate  u 
Lower  and  upper  limits  for  the  surface  coordinate  v 

Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 
region 

Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 
normal 

Components  of  the  translation  vector  from  the  global  origin  to  the  local 
origin  (default  is  ti  =  t2  =  t3  =  0) 


Surface  definition: 


The  local  coordinate  system  must  be  set  up  so  that  the  axis  of  the  cylinder  coincides 
with  the  2-axis.  In  this  case  u  —  z  and  v  is  related  to  the  angle  6  from  the  x-axis  in  the  xy 
plane  by 


v  =  tan  1 


( 


^  tan#^ 


The  local  coordinates  and  normal  vector  are  related  to  the  u,  v  coordinates  by 


x  =  a  cos  v 

— 

y  =  b  sinu 

ny  = 

£ 

II 

nz  = 

bcosv 

yfa?  —  (a2  —  62)  cos2  v 
a  sinn 

yfa?'  —  (a2  —  W)  cos2  v 
0 
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Equal  increments  in  v  are  not  equal  increments  in  9,  and  vise  versa.  However,  quadrants 
in  the  xy-plane  do  come  out  the  same  as  shown  in  the  following  table: 


V 

X 

y 

9 

0 

a 

0 

0 

7t/2 

0 

b 

7t/2 

7 r 

—a 

0 

7 r 

37t/2 

0 

-b 

37t/2 
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sphere 

spherical  surface 


Command:  sphere,  a 

Ul,  Uu,  Vl ,  Vu 
i,  Nvi  Qui  Qv 
h,  [^lj  t2,  <3] 


a 


Radius  of  sphere 


Ul,  Ujj  Lower  and  upper  limits  for  the  surface  coordinate  u 

Vl,  Vu  Lower  and  upper  limits  for  the  surface  coordinate  v 


N  N 

iyUl  iyV 


Qu,  Q 


V 


Iz 


t\,  h,  h 


Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 
region 

Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 
normal 

Components  of  the  translation  vector  from  the  global  origin  to  the  local 
origin  (default  is  ti  =  t2  =  t3  =  0) 


Surface  definition: 

The  local  coordinate  system  must  be  set  up  so  that  the  polar  axis  coincides  with  the 
axis.  In  this  case  u  is  the  angle  from  the  2-axis  and  v  is  the  angle  from  the  x-axis  in  the 
xy  plane.  The  local  coordinates  and  normal  vector  are  related  to  the  u,  v  coordinates  by 


x  =  a  sin  u  cos  v 

nx  =  sin  u  cos  v 

y  =  a  sin  u  sin  v 

ny  =  sin  u  sin?; 

z  =  a  cos  u 

nz  =  cos  u 
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oblatespheroid 

oblate  spheroidal  surface  obtained  by  rotating  ellipse  around  minor  axis 


Command:  oblate_spheroid,  a,  b 

Ul,  Uu,  Vl,  Vu 
Nui  Qi i)  Qv 

h,  [^i>  h,  ^3] 


a 

b 

uL,  uv 

vL,  vv 
Nu,  Nv 


Qu,Q 


v 


h 


^1)  ^2)  tz 


Half  the  major  diameter  of  elliptical  cross  section 
Half  the  minor  diameter  of  elliptical  cross  section 
Lower  and  upper  limits  for  the  surface  coordinate  u 
Lower  and  upper  limits  for  the  surface  coordinate  v 

Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 
region 

Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 
normal 

Components  of  the  translation  vector  from  the  global  origin  to  the  local 
origin  (default  is  ti  =  <2  =  £3  =  0) 


Surface  definition: 


The  local  coordinate  system  must  be  set  up  so  that  the  axis 
the  2-axis.  In  this  case  v  is  the  angle  from  the  a:- axis  in  the  xy 
angle  6  from  the  2-axis  by 

u  —  tan-1 


of  rotation  coincides  with 
plane,  u  is  related  to  the 
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The  local  coordinates  and  normal  vector  are  related  to  the  u,  v  coordinates  by 

x  =  asinucost; 

y  =  asintf  sini; 
z  =  b  cos  u 

Equal  increments  in  u  do  not  produce  equal  increments  in  6,  and  vise  versa.  However, 
quadrants  do  come  out  the  same  as  shown  in  the  following  table: 

'  u  x  y  z  6 

0  0  0  6  0 

7r/2  acosv  asinu  0  7r/2 

7T  0  0  — 6  7T 
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prolatespheroid 

prolate  spheroidal  surface  obtained  by  rotating  ellipse  around  major  axis 

Command:  prolate_spheroid,  a,  b 

Ul,  Uu,  Vl,  Vu 

Nu,  Qui  Qv 

h-,  [£i5  h,  h] 

a  Half  the  major  diameter  of  elliptical  cross  section 

b  Half  the  minor  diameter  of  elliptical  cross  section 

uL,  uu  Lower  and  upper  limits  for  the  surface  coordinate  u 

vL,  vu  Lower  and  upper  limits  for  the  surface  coordinate  v 

Nu,  Nv  Number  of  subdivisions  in  the  u  and  v  directions  in  the  planar  rectangular 

region 

Qui  Qv  Number  of  Gaussian  quadrature  points  per  subarea  in  the  u  and  v  directions 

Jz  Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 

normal 

t\,  t2,  £3  Components  of  the  translation  vector  from  the  global  origin  to  the  local 

origin  (default  is  £1  =  £2  =  £3  —  0) 

Surface  definition: 

The  local  coordinate  system  must  be  set  up  so  that  the  axis  of  rotation  coincides  with 
the  2- axis.  In  this  case  v  is  the  angle  from  the  x-axis  in  the  xy  plane  and  u  is  related  to  the 
angle  6  from  the  2-axis  by 
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The  local  coordinates  and  normal  vector  are  related  to  the  u,  v  coordinates  by 


,  .  a  sin  w  cosy 

x  =  b  sm  u  cos  v  nx  =  — — .  — 

-^/a2  —  (a2  —  b2)  cos2  v 
a  sin  u  sin  v 
2  —  (a2  —  b 2)  cos2  u 
bcosu 

z  =  a  cos  u  nz  =  —  — 

yjb2  —  (a2  —  W)  cos2  u 

Equal  increments  in  u  are  not  equal  increments  in  6,  and  vise  versa.  However,  quadrants 
do  come  out  the  same  as  shown  in  the  following  table: 

u  x  y  z  6 

0  0  0  a  0 

7t/2  b  cosy  6siny  0  7r/2 

7r  0  0  —a  7 r 


y  =  osinusiny 


ny  — 


yfb 
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fe_axisymmetric 

curve  generated  by  finite  element  interpolation  functions  is  rotated  around  z-axis 

Command:  fe_axisymmetric,  file.crd,  file.elt,  Qu,  Qv 

h ;  [ti,  h,  h] 

file.crd  File  containing  nodal  coordinates.  This  file  is  an  ASCII  file  containing  three 

columns  of  data  with  one  or  more  spaces  between  columns.  The  first  column 
contains  nodal  numbers,  the  second  column  contains  the  radial  r  coordi¬ 
nates  of  the  nodes,  and  the  third  column  eontans  the  axial  2  coordinates  of 
the  nodes.  The  nodes  lie  along  the  generating  curve  for  the  surface.  This 
curve  is  rotated  around  the  2- axis  to  generate  the  surface. 

file.elt  File  containing  element  node  numbers.  This  file  is  an  ASCII  file  that 

contains  two  columns  if  linear  interpolation  is  used  and  three  columns  if 
quadratic  interpolation  is  used.  The  columns  are  separated  by  one  or  more 
spaces.  Each  line  of  data  contains  the  node  numbers  for  an  element.  If 
there  are  two  columns,  the  node  numbers  are  for  the  first  and  last  node  in 
the  element.  If  there  are  three  columns,  then  the  columns  correspond  to 
first  node,  middle  node,  and  last  node  respectively. 

Qu,  Qv  Number  of  Gaussian  quadrature  points  in  the  u  and  v  directions 

Iz  Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 

normal 

t\,  ^2>  ts  Components  of  the  translation  vector  from  the  global  origin  to  the  local 

origin  (default  is  ti  =  t2  =  t3  =  0) 
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Surface  definition: 


The  local  coordinate  system  is  set  up  with  the  2-axis  being  the  axis  of  rotation.  Each 
finite  element  corresponds  to  one  CHIEF  subdivision  of  the  surface.  The  surface  coordinate 
v  is  the  angle  from  the  sr-axis  in  the  ary-plane.  The  surface  coordinate  u  is  in  the  interval 
[-1,1].  The  local  coordinates  and  normal  vector  are  related  to  the  u,  v  coordinates  by 


x  =  r(w)  cos  v 
y  =  r{u )  sinn 
2  =  z{u) 


nx  =  nr(u)  cosv 
ny  =  nr(u)  sinr; 
n2  =  nz(u) 


where 


and 


tangent  vector 


(nr(u) 
nz(u) 


-tz(u)/ 

tr(u)/^t2r(u)  +  t2z(u)  ) 


rn,  zn  given  coordinates  of  nodes  belonging  to  element 

4>n{u)  finite  element  interpolation  function  (linear  or  quadratic) 

For  linear  interpolation,  the  interpolation  functions  (j>n  are  given  by 

^i(w)  =  J(1  -«) 

Mu)  =  |(1  +«)• 

For  quadratic  interpolation, 

=  \u{u  -  1) 

(j>2  (u)  —  1  ~  U2 
4>z{u)  =  \u(u  -1-  1). 

The  normal  has  been  chosen  so  that  it  points  to  the  left  as  you  travel  along  the  cross- 
section  curve. 
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fe_3d 

surface  generated  by  finite  element  interpolation  functions 

Command:  fe_3d,  file.crd,  file.elt,  Qu,  Qv 

Iz,  [*i>  h,  £3] 

file.crd  File  containing  nodal  coordinates.  This  file  is  an  ASCII  file  containing  four 

columns  of  data  with  one  or  more  spaces  between  columns.  The  first  column 
contains  nodal  numbers  and  the  second  through  fourth  columns  contain  the 
x,  y,  z  coordinates  of  the  nodes. 

file.elt  File  containing  element  node  numbers.  This  file  is  an  ASCII  file  that 

contains  four  columns  if  linear  interpolation  is  used  and  eight  columns  if 
quadratic  interpolation  is  used.  The  columns  are  separated  by  one  or  more 
spaces.  Each  line  of  data  contains  the  node  numbers  for  an  element.  If  there 
are  four  columns,  the  node  numbers  are  for  the  corner  nodes  of  the  element 
taken  in  a  counterclockwise  direction  around  the  element.  If  there  are  eight 
columns,  then  the  first  four  columns  are  the  same  as  for  the  linear  case 
and  the  last  four  columns  are  for  the  midnodes  taken  in  a  counterclockwise 
direction  beginning  with  the  node  between  the  first  and  second  nodes. 

Qv  Number  of  Gaussian  quadrature  points  in  the  u  and  v  directions 

Iz  Specifies  rotation  from  local  to  global  coordinate  system  and  orientation  of 

normal 

tu  t2,  £3  Components  of  the  translation  vector  from  the  global  origin  to  the  local 

origin  (default  is  £1  =  £2  =  £3  =  0) 
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Surface  definition: 


Each  finite  element  corresponds  to  one  CHIEF  subdivision  of  the  surface.  Figure  5  shows 
the  base  rectangle  in  the  uv  coordinate  system. 


V 
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Figure  5.  The  base  rectangle  in  the  uv  coordinate  system. 


The  local  coordinates  and  normal  vector  are  related  to  the  u,  v  coordinates  by 


x  =  x(u,  v ) 
V  =  y(«>  v) 
z  =  z(u,  V ) 

where 


nx  =  tx(u,  v)/ tju,  u)2  +  tju,  v )2  +  tju,  v )2 
ny  =  ty(u,  v)/  J tju,  v)2  +  ~tju,  u)2  +  ~tju,  u)2 
nz  =  tju,  v)/ tju,  v)2  +  ty(u,  v)2  +  tju,  v)2, 


^ x(u ,  v)^ 

4,8 

(x\ 

y(u,  v ) 

=  J24>n(u,v) 

Vn 

n~l 

\Zn) 
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and 


v)^ 

Y^d(f)n{u,v) 

du 

- 1 

fx\ 

ty(u,  V ) 

Vn 

v) ) 

71 — 1 

\Zn) 

vj\ 

tVy(U,  V ) 

d(f>n(u,  v ) 

2-*  dv 

(x\ 
•*- n 

Vn 

\tVz(u>v)) 

71=1 

\Zn) 

u  tangent  vector 


v  tangent  vector 


tx  (u, 

ty(u,V) 

\tz(u,v)j 


(«,  ^(ti,  «)  -  «)tJ(U,  V)^ 

*“(«,  «)*£(«>  V ) 
^(U,  V)tVy(U,  V )  -  tUy(U,  V)tl(U,  V )  } 


Xu-,  Vni  zn  given  coordinates  of  nodes  belonging  to  element 
(pn(u,v)  finite  element  interpolation  function  (linear  or  quadratic) 

For  linear  interpolation,  the  interpolation  functions  0n  are  given  by 

<t>\{u,v)  =  i(l-u)(l  +  u) 

<t>2{u,v)  =  |(1  ~u)(l  -v) 

<t>z{u,v)  =  j(l  +  w)(l  —  v) 

<j>4(u,v)  =  j(l  +  «)(1  +  V). 

For  quadratic  interpolation, 

=  |(1  -  w)(l  +  v)(u-w-  1) 

(j)2(u, v )  =  \(1  -u)(  1  -  v)(-v  -u-  1) 

<p3(u,  v)  =  |(1  +  u)(  1  -v)(-v  +  u-  1) 

04 (u,  v )  =  |(1  +  u)(l  +  i;)(u  +  u  -  1) 

05 («,«)  =  |(1  -M)(l+U)(l  -V) 

<j)6{u,  v)  =  |(1  -  u)(l  -  v)(l  +  u ) 

4>7(u,  v)  =  |(1  +  «)(1  -  v)(l  +  v ) 
c/)8(u,  v )  =  |(1  +  u)(l  +  t^)(l  -  u ). 


The  normal  was  chosen  so  that  the  nodes  appear  in  counterclockwise  order  when  viewed 
from  the  positive  normal  side  of  the  surface. 
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5.  CONTROL  SECTION 


In  this  section,  commands  are  given  to  specify  the  type  of  job  CHIEF  is  asked  to  perform. 
The  commands  in  this  section  may  be  given  in  any  order.  Only  the  first  three  characters  of 
the  command  word  itself  are  used  by  the  program.  Other  forms  of  the  command  words  are 
perfectly  valid  as  long  as  the  same  first  three  characters  are  used. 

Each  command  line  consists  of  a  command  or  a  command  followed  by  one  or  more 
arguments.  Some  command  lines  must  be  followed  by  additional  lines  of  data.  Commands, 
arguments,  and  data  items  are  separated  by  commas.  Parameters  that  have  been  defined 
in  the  parameter  section  and  expressions  involving  these  parameters  may  be  used  for  any 
numerical  input  in  this  section.  As  always,  blank  lines,  comments,  and  spaces  are  ignored. 

The  following  types  of  problems  may  be  run  in  CHIEF  (with  only  one  type  allowed  per 
run).  The  first  two  types  are  the  most  commonly  used.  Following  each  type  are  a  list  of 
commands  that  must  be  included  in  the  input  file.  Optional  commands  are  given  in  braces. 

1.  Calculate  surface  and/or  field  pressures  (see  example  problems  1,  2,  5,  6,  8): 

velocity  or  planewave 
surface  pressures 
[fax-field] 

[near-field] 

2.  Calculate  impedance  and/or  field  matrices  (see  example  problem  4): 

impedance 
[far-f  ield] 

[near-field] 

3.  Calculate  A  and  B  matrices  (see  example  problem  9): 

a_b  matrices 

One  or  more  of  the  following  are  required: 

[surface  pressures] 

[fax-field] 

[neax-field] 

4.  Point  source  check  on  numerical  accuracy  (see  example  problem  3): 

point  souxce 
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CHIEF  solves  exterior  radiation  or  scattering  problems  by  default.  If  the  inside  com¬ 
mand  is  given,  CHIEF  solves  an  interior  problem.  For  interior  problems,  commands  such 
as  far-field  and  infinite  baffle  make  no  sense.  We  now  provide  the  format  for  each 
command.  Optional  parts  of  the  command  line  are  in  brackets. 
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a_b  matrices 

outputs  the  surface  or  field  A  and  B  matrices 


Command:  a_b  matrices 

When  the  surface  pressures  command  is  given,  CHIEF  produces  a  system  of  equations 
of  the  form  AP  =  BV  +  Pinc,  where  P  is  the  vector  of  surface  pressures,  V  is  the  vector  of 
surface  normal  velocities,  and  Pinc  is  the  vector  of  incident  surface  pressures  resulting  from 
an  incident  plane  wave  specified  by  the  user.  When  the  fax-field  pressures  command  is 
given,  CHIEF  produces  a  set  of  equations  of  the  form  Pff  =  AP  +  BV ,  where  P//  is  the 
vector  of  far-field  radiated  pressures  at  prescribed  angles.  When  the  near-field  pressures 
command  is  given,  CHIEF  produces  a  set  of  equations  of  the  form  Pnf  =  AP  +  BV ,  where 
Pnf  is  the  vector  of  near-field  pressures  at  prescribed  locations. 

A  and  B  matrices  may  be  output  for  any  or  all  of  the  above  three  cases.  These  matrices 
are  output  to  binary  files  with  extensions  .  abs,  .  abf ,  and/or  .  abn  if  the  commands  surface 
pressures,  far-field  pressures,  and/or  near-field  pressures,  respectively,  are  in  the 
CHIEF  input  file. 
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far-field 

calculates  far-field  pressures  or  matrix 


Command: 


far-field,  N,  Iprt 

Q\nci  &L1  &U )  tfiinc 


(in  ai\  nj\  x/v  J./V 

°L  )  Vui  °ina  <Pl>  %>  <Pinc 


N 


Number  of  far-field  data  lines  to  follow  command 


Iprt 


Ov,  Oinc 

^Li  4>Ui  fiinc 


Print  flag.  If  Iprt=l,  then  far-field  pressures  are  written  to  the  standard 
output  file  (< input  f  ilename> .  out).  If  Iprt=2,  then  the  far-field  pressures 
are  written  to  an  ASCII  file  (<input  filename>.ffp)  containing  only  far- 
field  output.  If  Iprt=3,  then  the  far-field  pressures  are  written  to  both  the 
standard  output  file  and  the  far-field  ASCII  output  file.  If  far-field  matrices 
are  computed,  they  are  written  to  the  binary  file  (cinput  f  ilename>.pmf) 
regardless  of  the  value  of  the  print  flag. 

Lower  limit,  upper  limit,  and  increment  in  0  for  k- th  0,  <j)  range 
Lower  limit,  upper  limit,  and  increment  in  4>  for  k- th  0,  4>  range 


Each  data  line  specifies  a  set  of  equally  spaced  0’s  and  a  set  of  equally  spaced  0’s.  The 
program  performs  far-field  computations  for  all  0,  4>  pairs,  where  6  is  taken  from  the  set  of 
equally  spaced  0’s  and  <p  is  taken  from  the  set  of  equally  spaced  4>'s.  Each  successive  data  line 
specifies  another  collection  of  0,  (f>  pairs  for  which  far-field  calculations  are  to  be  performed. 

If  the  Impedance  command  is  contained  in  the  control  section,  then  fax-field  matrices  axe 
computed.  Each  matrix  row  corresponds  to  a  direction  specified  by  a  0,  0  pair,  and  each 
column  corresponds  to  a  surface  subdivision.  If  G  is  the  fax-field  matrix,  then  the  vector  of 
far-field  pressures  Pjj  at  the  specified  directions  is  related  to  the  vector  of  surface  normal 
velocities  V  by 

Pff  =  GV. 

If  a  surface  velocity  distribution  is  specified  by  a  Velocity  command  in  the  control  section, 
then  the  far-field  radiated  pressures  corresponding  to  the  velocity  distribution  are  computed 
instead  of  the  far-field  matrix.  If  an  incident  plane  wave  is  specified  through  the  Plane  wave 
command,  then  the  scattered  far-field  pressures  are  computed. 
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fluid 


user  specified  fluid  density  and  sound  speed 


Command:  fluid,  p,  c 

p  Density  of  fluid 

c  Sound  speed  in  fluid 

If  the  fluid  command  is  omitted,  then  the  default  values  of  p  =  1000  and  c  =  1500  will 

be  used.  These  values  are  fairly  typical  for  water  in  mks  units,  but  if  a  different  fluid  or  a 
different  set  of  units  is  used,  then  the  appropriate  values  can  be  set  with  this  command. 
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frequency 

specifies  frequency  range  for  calculations 
Command:  frequency,  fL  [,  fv,  finc] 

fh  Lower  frequency  limit  for  range  of  equally  spaced  frequencies 

fu  Upper  frequency  limit  for  range  of  equally  spaced  frequencies 

fine  Frequency  increment 

If  fu  and  finc  are  missing,  then  only  one  frequency  /L  will  be  used. 
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impedance 

calculates  impedance  matrix 


Command:  impedance,  Isym,  Iprt 


Isym 


Iprt 


If  Isym  =  1,  then  the  velocity  distribution  has  the  same  symmetry  as  the 
radiating  surface.  If  Isym  =  0,  then  CHIEF  will  not  take  advantage  of 
velocity  symmetry 

Print  flag.  If  Iprt  =  1,  then  the  impedance  matrix  will  be  output  to  the 
binary  file  < input  filename>.zmx. 


The  impedance  matrix  Z  is  defined  by 


F  =  ZV, 


where  V  is  the  vector  of  surface  normal  velocities  on  the  CHIEF  subdivisions  and  F  is  the 
vector  of  surface  forces  on  the  CHIEF  subdivisions. 
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infinite  baffle 

specifies  infinite  rigid  baffle 


Command:  infinite  baffle 

CHIEF  treats  the  radiating  surface  as  a  planar  region  in  an  infinite  rigid  baffle.  The  surface 
definition  must  only  consist  of  planar  regions. 
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inside 

specifies  that  an  interior  problem  is  to  be  solved 


Command:  inside 

CHIEF  solves  an  interior  problem  instead  of  an  external  problem.  The  positive  normal 
direction  is  still  taken  as  external  to  the  surface. 
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interior  points 

adds  equations  corresponding  to  prescribed  interior  points  to  improve  numerical  stability  of 
CHIEF  equations  at  certain  critical  frequencies 


Command:  interior  points,  N 

*1,  2/i,  zi 

XN,  Vn,  Zn 

N  Number  of  interior  points  per  symmetry  block 

xki  Vk,  zk  Global  cartesian  coordinates  of  k- th  interior  point  where  k  =  1,...,  AT. 

The  CHIEF  surface  A  matrix  becomes  ill-conditioned  near  certain  critical  frequencies 
corresponding  to  resonances  of  the  interior  problem  with  a  zero  surface  pressure  bound¬ 
ary  condition.  The  addition  of  equations  corresponding  to  interior  points  relieves  this  ill- 
conditioning  if  at  least  one  of  the  points  does  not  lie  on  a  nodal  surface  for  the  interior  mode 
corresponding  to  the  critical  frequency. 


44 


near-field 

calculates  near-field  pressures  or  matrix 


Command:  near-field,  N,  Iprt 

*i»  2/i)  *i 


VNi  zn 


N 

Iprt 


Xk,  Vk,  zk 


Number  of  near-field  points 

Print  flag.  If  Iprt  =  1,  then  near-field  pressures  are  written  to  the  stan¬ 
dard  output  file  (<input  filename>.out).  If  Iprt  =  2,  then  the  near-field 
pressures  are  written  to  an  ASCII  file  (<input  f ilename> .nfp)  contain¬ 
ing  only  near-field  output.  If  Iprt  =  3,  then  the  near-field  pressures  are 
written  to  both  the  standard  output  file  and  the  near-field  ASCII  output 
file.  If  near-field  matrices  are  computed,  they  are  output  to  the  binary  file 
(<input  f  ilename>  .pmn)  regardless  of  the  value  of  the  print  flag. 

Global  cartesian  coordinates  of  fc-th  near-field  point  at  which  near-field 
pressures  are  desired  (k  =  l, ...  ,N). 


If  the  Impedance  command  is  contained  in  the  control  section  ,  then  near-field  matrices 
are  computed.  Each  row  of  this  matrix  corresponds  to  a  near-field  point  and  each  column 
corresponds  to  a  surface  subdivision.  If  G  is  the  near-field  matrix,  then  the  vector  of  near¬ 
field  pressures  Pnf  at  the  specified  points  is  related  to  the  vector  of  surface  normal  velocities 
V  by 

Pnf  =  GV. 

If  a  surface  velocity  distribution  is  specified  by  a  Velocity  command  in  the  control  section, 
then  the  near-field  radiated  pressures  corresponding  to  the  velocity  distribution  are  computed 
instead  of  the  far-field  matrix.  If  an  incident  plane  wave  is  specified  through  the  Plane  wave 
command,  then  the  scattered  near-field  pressures  are  computed. 
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plane  wave 

specifies  rigid  scattering  problem 


Command:  plane  wave,  Isym,  Pinc,  6inc,  <j>inc 


Isym 


p. 

1  me 

@inc  ■ 


0inc 


If  Isym  —  1,  then  the  incident  pressure  field  has  the  same  symmetry  as 
the  radiating  surface.  If  Isym  =  0,  then  CHIEF  will  not  take  advantage  of 
velocity  symmetry 

Amplitude  of  the  incident  plane  wave 

Angles  (in  degrees)  specifying  the  direction  from  which  the  plane  wave  is 
coming  relative  to  the  global  coordinate  axes.  9inc  is  the  angle  from  the 
rE3-axis  and  <pinc  is  the  angle  from  the  Xi-axis  in  the  ^i^-plane. 


Since  CHIEF  does  not  contain  any  structural  capabilities,  it  can  only  compute  rigid  scatter¬ 
ing. 
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point  source 

uses  pressures  from  specified  point  sources  to  check  the  accuracy  of  CHIEF  calculations 


Command: 


point  source,  Isym,  N 

*1,  Vu  Zi,  wr,  A“‘ 


xN,  V„,  is,  »'v.  »T.  if 

Isym  =  1  specifies  that  the  point  source  distribution  has  the  same  symmetry 
as  the  radiating  surface.  In  this  case,  the  point  sources  only  need  to  be  input 
for  the  first  symmetry  block 

Number  of  point  sources 

Global  cartesian  coordinates  of  k- th  point  source 

Real  and  imaginary  parts  of  the  weight  for  the  fc-th  point  source 

Flag  that  indicates  whether  the  k- th  point  source  is  interior  to  (I%.xt  =  0) 
or  exterior  to  {I%xt  =  1)  the  radiating  surface 

The  point  source  check  uses  analytic  expressions  to  calculate  surface  velocities  and  in¬ 
cident  surface  pressures  (for  exterior  point  sources)  due  to  the  prescribed  point  sources 
assuming  that  the  surface  is  an  imaginary  surface  in  the  free  field.  The  surface  pressures  are 
then  calculated  using  the  CHIEF  equations 

AP  =  BV  +  Pinc. 

The  surface  pressures  calculated  using  CHIEF  can  be  compared  to  the  pressures  on  the 
surface  calculated  using  analytical  expressions  for  point  sources.  This  comparison  can  be 
used  to  determine  the  adequacy  of  the  surface  break-up  into  subareas,  adequacy  of  interior 
points  to  remove  ill-conditioning  at  critical  frequencies,  and  errors  due  to  improper  input  such 
as  wrong  directions  for  normal  vectors.  The  location  of  the  point  sources  should  be  selected 
to  produce  a  velocity  distribution  that  is  as  smooth  as  the  expected  velocity  distribution. 
For  thin  bodies,  the  velocity  due  to  an  interior  point  source  will  vary  very  rapidly  near  the 
point  source.  Therefore,  it  may  be  necessary  to  use  exterior  point  sources  to  obtain  a  velocity 
distribution  with  the  desired  smoothness. 


Isym 


N 

xk,  Vk,  *k 

wr,  wim 

jext 
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rhs  symmetry 

specifies  that  the  velocity  (appearing  on  the  right-hand  side  of  the  equation  Ap  =  Bv)  has 
the  same  symmetry  as  the  surface 


Command:  rhs 
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surface  pressures 

specifies  calculation  of  surface  pressures 


Command: 

Iprt 


surface  pressures,  [Iprt] 

Print  flag.  If  Iprt  =  1,  then  surface  pressures  are  written  to  the  standard 
output  file  (Cinput  f  ilename> .  out).  If  Iprt  =  2,  then  the  surface  pressures 
are  written  to  an  ASCII  file  (cinput  filenames  sps)  that  only  contains 
surface  pressure  output.  If  Iprt  =  3,  then  the  surface  pressures  are  written 
to  both  the  standard  output  file  and  the  surface  pressure  ASCII  output  file. 
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velocity 

specifies  surface  normal  velocity  input 


Command:  velocity,  <filename> 

<f  ilename>  Name  of  file  containing  surface  normal  velocity  input.  This  file  is  an  ASCII 
file  containing  five  columns  separated  by  one  or  more  spaces.  The  first  col¬ 
umn  is  the  velocity  distribution  number  (i.e.,  the  first  velocity  distribution 
will  have  a  1  in  this  column,  the  second  a  2,  etc.).  The  second  column 
contains  the  symmetry  block  number.  The  third  column  contains  the  sur¬ 
face  subdivision  number.  The  fourth  and  fifth  columns  contain  the  real  and 
imaginary  parts  of  the  normal  velocity. 

If  the  velocity  has  the  same  symmetry  as  the  radiating  surface,  then  only  one  symmetry 
block  of  the  velocities  needs  to  be  input. 
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velocity  (alternate  form) 

specifies  surface  normal  velocity  input  by  means  of  a  loop 


Command:  vel(ii  :  i2  :  *3,  ji  :  j2  :  j3,  kx  :  k2  :  A:3)  =  ( x ,  y)  I  —  h,  I2,  h 


I 


h,  h,  h 


i\  :  h  ■  h 


j  1  ;  32  '■  jz 


ki  :  k2  :  k3 


x,y 


loop  index  (does  not  have  to  be  I — the  name  can  be  selected  by  the  user). 
If  no  loop  is  desired  the  portion  of  the  command  beginning  with  “J  =”  can 
be  omitted. 

Ii  is  the  starting  index  for  the  loop,  I2  is  the  ending  index  for  the  loop,  and 
I3  is  the  step  size.  These  limits  can  be  replaced  by  any  expression  involving 
previously  defined  parameters. 

specifies  a  range  for  the  surface  subdivision  number,  where  i\  is  the  starting 
index,  i2  is  the  ending  index,  and  i3  is  the  step  size.  If  iz  is  not  present, 
it  is  set  to  1.  If  i2  and  i3  are  not  present,  i2  is  set  equal  to  ix.  The  colon 
is  not  necessary  if  the  following  index  is  not  present  (e.g.,  the  range  could 
be  given  by  i\  :  i2  or  just  ix).  The  range  limits  can  be  replaced  by  any 
expression  involving  previously  defined  parameters  and/or  the  index  I. 

specifies  a  range  for  the  symmetry  block  number,  where  j\  is  the  starting 
index,  j2  is  the  ending  index,  and  jz  is  the  step  size.  If  j3  is  not  present, 
it  is  set  to  1.  If  j2  and  j3  are  not  present,  j2  is  set  equal  to  ji-The  colon 
is  not  necessary  if  the  following  index  is  not  present  (e.g.,  the  range  could 
be  given  by  j\  :  j2  or  just  ji).  The  range  limits  can  be  replaced  by  any 
expression  involving  previously  defined  parameters  and/or  the  index  I. 

specifies  a  range  for  the  velocity  distribution  number,  where  Aq  is  the  start¬ 
ing  index,  k2  is  the  ending  index,  and  k3  is  the  step  size.  If  k3  is  not  present, 
it  is  set  to  1.  If  k2  and  k3  are  not  present,  k2  is  set  equal  to  fci.The  colon 
is  not  necessary  if  the  following  index  is  not  present  (e.g.,  the  range  could 
be  given  by  ki  :  k2  or  just  kx).  The  range  limits  can  be  replaced  by  any 
expression  involving  previously  defined  parameters  and/or  the  index  I. 

x  is  the  real  part  of  the  velocity  and  y  is  the  imaginary  part  of  the  velocity. 
If  the  imaginary  part  is  zero,  the  user  can  use  x  in  place  of  (x,  y).  x  and  y 
can  be  replaced  by  any  expression  involving  previously  defined  parameters 
and/or  the  index  I. 


If  the  velocity  has  the  same  symmetry  as  the  radiating  surface,  then  only  one  symmetry 
block  of  the  velocities  needs  to  be  input. 
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6.  OUTPUT  FILE  FORMATS 


In  this  section,  we  describe  the  formats  of  the  various  output  files.  The  files  that  are 
likely  to  be  small  are  written  as  formatted  ASCII  files.  The  files  that  are  likely  to  be  large  are 
written  as  unformatted  binary  files.  The  version  of  this  program  for  Windows  95/98®  and 
Windows  NT®  machines  was  compiled  using  the  Lahey  Fortran  95  compiler.  The  version 
for  the  DEC  Alpha  OSF1  was  compiled  using  the  Digital  Fortran  90  compiler.  In  both  cases, 
the  unformatted  sequential  binary  files  are  divided  into  records  that  have  a  4-byte  integer 
prefix  and  suffix  containing  the  number  of  bytes  in  the  record  (excluding  prefix  and  suffix) . 
Both  also  use  IEEE  little  endian  format  for  floating  point  reals  and  integers.  The  output  at 
each  frequency  in  a  binary  file  begins  with  a  standard  header  record  containing  the  following 
12  items: 

1.  Frequency 

2.  Fluid  density 

3.  Fluid  sound  speed 

4.  Symmetry  type  (’ROT’,  ’REF’,  or  ’NONE’) 

5.  Number  of  symmetry  blocks 

6.  Number,  of  symmetry  blocks  output 

7.  Velocity  symmetry  flag(0  =  no  symmetry,  1  =  same  symmetry  as  surface) 

8.  Number  of  rows  in  A  and  B  matrices  (number  of  areas  plus  number  of  interior  points) 

9.  Number  of  areas  (i.e.,  total  number  of  subdivisions  for  radiating  surface) 

10.  Number  of  theta-phi  field  point  pairs 

11.  Number  of  near-field  points 

12.  Number  of  velocity  distributions. 

The  formats  of  quantities  in  the  various  output  files  are  described  using  Fortran  type  speci¬ 
fiers.  For  example,  (F12.3)  specifies  a  fixed  point  style  real  number  with  a  total  width  of  12 
and  3  digits  to  the  right  of  the  decimal  point,  (15)  specifies  a  five  digit  integer,  and  (E14.2) 
specifies  an  exponential  style  real  number  with  a  total  width  of  14  and  3  digits  to  the  right  of 
the  decimal  place.  The  total  width  includes  any  preliminary  plus  or  minus  sign,  the  decimal 
point,  and  the  exponent  (E±xx).  We  will  now  look  at  each  file  type  individually. 
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Surface  pressure  file  (.sps) 

This  is  a  formatted  ASCII  file.  The  first  line  contains  the  frequency  with  the  format 
(F12.3).  The  following  lines  have  the  form, 


irhs  hik  i  Preai{i)  Pimag(i )  area(i), 


where  irfls  is  an  index  running  over  the  number  of  velocity  distributions,  i&zfc  is  an  index 
running  over  the  number  of  symmetry  blocks,  i  is  an  index  running  over  the  number  of 
CHIEF  subdivisions,  PTeai{i)  and  Pimag  are  the  real  and  imaginary  parts  of  the  surface 
pressure  on  the  i-th  subdivision,  and  area(i)  is  the  area  of  the  i-th  subdivision.  The  format 
for  each  of  these  lines  is  (315,  3E15.6).  The  index  i  varies  the  fastest,  next  is  IUk,  and  slowest 
is  Irhs.  The  above  lines  are  repeated  for  each  frequency. 

Far-field  pressure  file  (.ffp) 

This  is  a  formatted  ASCII  file.  The  first  line  contains  the  frequency  with  the  format 
(F12.3).  The  following  lines  have  the  form, 


irhs  <t>  e  Prff(<j>,9)  P}™9(4>,  9), 


where  irhs  is  an  index  running  oyer  the  number  of  velocity  distributions,  <p  and  9  run  over 
all  9-<t>  pairs,  and  Pjefal(4>,  9),  P™a9((f>,9)  are  the  real  and  imaginary  parts  of  the  far-field 
pressure.  The  format  for  each  of  these  lines  is  (15,  2F10.2,  2E15.6).  The  above  lines  are 
repeated  for  each  frequency. 

Near-field  pressure  file  (.nfp) 

This  is  a  formatted  ASCII  file.  The  first  line  contains  the  frequency  with  the  format 
(F12.3).  The  following  lines  have  the  form, 


*  y  z  P™fal{x,  y,  z)  P™9(x,y,z), 

where  irhs  is  an  index  running  over  the  number  of  velocity  distributions,  x,  y,  and  z  run  over 
all  prescribed  near-field  points,  and  P%?l(x,y,  z),  P™a9(x,y,z )  are  the  real  and  imaginary 
parts  of  the  near-field  pressure.  The  format  for  each  of  these  lines  is  (15,  3E12.4,  2E15.6). 
The  above  lines  are  repeated  for  each  frequency. 

Impedance  Matrix  file  (.zmx) 

This  is  a  binary  file.  The  first  record  is  the  standard  header  record  for  the  first  frequency. 
The  second  record  contains  the  surface  areas  of  the  CHIEF  subdivisions.  The  next  records 
contain  the  necessary  symmetry  blocks  of  the  impedance  matrix  (one  block  if  there  is  velocity 
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symmetry).  For  each  succesive  frequency,  there  is  a  header  record  followed  by  the  necessary 
blocks  of  the  impedance  matrix.  The  matrices  are  written  out  column  by  column. 

Far-field  Matrix  file  (.pmf) 

This  is  a  binary  file.  The  first  record  is  the  standard  header  record  for  the  first  frequency. 
The  second  record  contains  the  theta-phi  pairs  for  the  far-field  directions.  The  next  records 
contain  the  necessary  symmetry  blocks  of  the  far-field  matrix  (one  block  if  there  is  velocity 
symmetry).  For  each  succesive  frequency,  there  is  a  header  record  followed  by  the  necessary 
blocks  of  the  far-field  matrix.  The  matrices  are  written  out  column  by  column. 

Near-field  Matrix  file  (.pmn) 

This  is  a  binary  file.  The  first  record  is  the  standard  header  record  for  the  first  frequency. 
The  second  record  contains  the  coordinates  for  the  near-field  points.  The  next  records 
contain  the  necessary  symmetry  blocks  of  the  near-field  matrix  (one  block  if  there  is  velocity 
symmetry) .  For  each  succesive  frequency,  there  is  a  header  record  followed  by  the  necessary 
blocks  of  the  near-field  matrix.  The  matrices  are  written  out  column  by  column. 

Surface  A,B  Matrices  file  (.abs) 

This  is  a  binary  file.  The  first  record  is  the  standard  header  record  for  the  first  frequency. 
The  second  record  contains  the  surface  areas  of  the  CHIEF  subdivisions.  The  next  records 
contain  the  necessary  symmetry  blocks  of  the  A,  B  matrices  (one  block  if  there  is  velocity 
symmetry).  For  each  symmetry  block,  the  A  matrix  is  output  first  to  a  record  and  the  B 
matrix  is  output  to  the  following  record.  For  each  succesive  frequency,  there  is  a  header 
record  followed  by  the  necessary  blocks  of  the  A,  B  matrices.  The  matrices  are  written  out 
column  by  column. 

Far-field  A,B  Matrices  file  (.abf) 

This  is  a  binary  file.  The  first  record  is  the  standard  header  record  for  the  first  frequency. 
The  second  record  contains  the  theta-phi  pairs  for  the  far-field  directions.  The  next  records 
contain  the  necessary  symmetry  blocks  of  the  A,  B  matrices  (one  block  if  there  is  velocity 
symmetry).  For  each  symmetry  block,  the  A  matrix  is  output  first  to  a  record  and  the  B 
matrix  is  output  to  the  following  record.  For  each  succesive  frequency,  there  is  a  header 
record  followed  by  the  necessary  blocks  of  the  A,  B  matrices.  The  matrices  are  written  out 
column  by  column. 

Near-field  A,B  Matrices  file  (.abn) 

This  is  a  binary  file.  The  first  record  is  the  standard  header  record  for  the  first  frequency. 
The  second  record  contains  the  coordinates  for  the  near-field  points.  The  next  records 
contain  the  necessary  symmetry  blocks  of  the  A,  B  matrices  (one  block  if  there  is  velocity 
symmetry) .  For  each  symmetry  block,  the  A  matrix  is  output  first  to  a  record  and  the  B 
matrix  is  output  to  the  following  record.  For  each  succesive  frequency,  there  is  a  header 
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record  followed  by  the  necessary  blocks  of  the  A,  B  matrices.  The  matrices  are  written  out 
column  by  column. 
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7.  EXAMPLE  PROBLEMS 


In  this  section,  we  will  show  how  to  set  up  the  ASCII  input  files  for  eight  example 
problems.  Most  of  these  examples  are  the  same  as  those  used  in  the  original  CHIEF  manual. 
To  gain  maximum  benefit  from  these  example  problems,  the  user  should  try  to  set  up  each 
of  the  input  files  without  looking  at  the  answers  provided. 
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Example  Problem  1 

Figure  6  shows  the  geometry  for  the  first  problem.  It  consists  of  a  free-flooded  cylinder 
having  an  inner  radius  of  1.68  m,  an  outer  radius  of  2.04  m,  and  a  height  of  1.2  m.  The 
normal  velocity  on  the  inner  cylindrical  surface  is  -1.029  m/sec,  the  normal  velocity  on  the 
outer  cylindrical  surface  is  0.971  m/sec,  and  the  normal  velocity  on  the  top  and  bottom 
surfaces  is  0.097  m/sec. 


x3 


9-12  are  on  the  inside  surface 


Figure  6.  Geometry  of  example  problem  1. 

The  problem  is  to  compute  the  surface  pressures  and  the  far-field  pattern  in  the  plane 
4>  =  0  with  Ad  =  5°  and  0  <  6  <  90°  for  a  frequency  of  1500/(27t).  Use  12-fold  rotational 
symmetry  and  no  velocity  symmetry  (even  though  there  really  is)  to  solve  the  problem.  The 
surface  and  far-field  pressures  should  be  output  to  a  file.  The  following  is  an  input  file  for 
this  problem: 
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freql=1500/(2*pi) 

nblks=12 

ai=l .68 

ao=2 . 04 

ht=l . 2 

htd2= . 5*ht 

rotlim=pi/nblks 

rho=999 . 89 

end*params 

rotational ,nblks 

circ_planax,  htd2 
ai , ao , -rotlim , rotlim 

2.1.4.16 
3, 0,0,0 

circ_cylinder ,  ao 

-htd2 , htd2 , -rotlim , rotlim 

4.1.4.16 
3 

circ_planar,  -htd2 
ai , ao , -rot lim, rotlim 

2.1.4.16 
-3 

circ_cylinder ,  ai 

-htd2 , htd2 , -rotlim , rotlim 

4.1.4.16 
-3 

end*geom 

fluid , rho , 1500 . 05 

freq,  freql,  freql,  0. 

vel,  tcl.vel 

surface  pressures, 3 

farfield  pressures, 1 ,3 
0,90,5,0,0,0 


The  user  also  needs  to  construct  a  file  containing  the  velocity  distribution.  Since  we  are 
not  using  the  symmetry  of  the  velocity,  the  velocity  distribution  must  be  input  for  all  12  of 
the  rotational  symmetry  blocks.  If  we  had  chosen  to  use  the  velocity  symmetry,  then  the 
input  file  would  only  need  to  contain  the  velocity  distribution  for  the  first  symmetry  block. 
The  velocity  input  file  for  example  problem  1  is  shown  below: 
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1 

4 
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1 

7  12 

-1.029 

0 

1 

8  1 

0.097 

0 

1 

8  2 

0.097 

0 

1 

8  3 

0.971 

0 

1 

8  4 

0.971 

0 

1 

8  5 

0.971 

0 

1 

8  6 

0.971 

0 

1 

8  7 

0.097 

0 

1 

8  8 

0.097 

0 

1 

8  9 

-1.029 

0 

1 

8  10 

-1.029 

0 

1 

8  11 

-1.029 

0 

1 

8  12 

-1.029 

0 

1 

9  1 

0.097 

0 

1 

9  2 

0.097 

0 

1 

9  3 

0.971 

0 

1 

9  4 

0.971 

0 

1 

9  5 

0.971 

0 

1 

9  6 

0.971 

0 

1 
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0 
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0 

1 
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0 

1 
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0 

1 

11  8 

0.097 

0 
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1 

11 

9 
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11 
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11 
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11 
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12 
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0.097 
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12 
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0.971 

0 

1 

12 
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0.971 

0 

1 

12 

5 

0.971 

0 

1 

12 

6 

0.971 

0 

1 

.12 

7 

0.097 

0 

1 

12 

8 

0.097 

0 

1 

12 

9 

-1.029 

0 

1 

12 

10 

-1.029 

0 

1 

12 

11 

-1.029 

0 

1 

12 

12 

-1.029 

0 

For  this  example,  the  velocity  could  be  specified  directly  in  the  input  file  instead  of 
reading  the  velocities  from  a  data  file.  The  line, 

vel,  tcl.vel 

should  be  replaced  by  the  lines, 

vel(l:2, 1:12, 1)=(. 097,0.) 
vel(3:6, 1:12, 1)=(. 971,0.) 
vel (7:8, 1:12, 1)=(. 097,0.) 
vel(9:12,l:12,l)=(-l. 029,0.) 

The  following  pages  show  the  output  that  would  be  obtained  in  the  .out  file  for  this 
problem. 
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Example  Problem  2 

Figure  7  shows  the  geometry  for  the  second  problem.  The  geometry  of  this  problem  is  the 
same  as  example  problem  1.  In  this  example,  three  planes  of  reflective  symmetry  are  used 
in  place  of  rotational  symmetry  and  velocity  symmetry  is  used.  The  surface  and  far-field 
pressures  should  be  output  to  a  file. 


Figure  7.  Geometry  of  example  problem  2. 
The  following  is  an  input  file  for  this  problem: 

f reql=1500/ (2*pi) 

ai=l .68 

ao=2 . 04 

ht=l .2 

htd2=.5*ht 

ref lim=pi/2 

end*params 

ref lective,3 

circ_planar,  htd2 
ai,ao,0,reflim 
2,3,4,16 
3, 0,0,0 

circ_cylinder,  ao 
0,htd2,0,reflim 
2,3,4,16 
3 
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circ_cylinder,  ai 
0,htd2,0,reflim 
2,3,4,16 
-3 

end*geom 

freq,  freql,  freql,  0. 

rhs  symmetry 

vel,  tc2.vel 

surface  pressures, 3 

farfield  pressures, 1,3 
0,90,5,0,0,0 


The  user  also  needs  to  construct  a  file  containing  the  velocity  distribution.  Since  we  chose 
to  use  the  velocity  symmetry,  the  input  file  only  needs  to  contain  the  velocity  distribution 
for  the  first  octant.  The  velocity  input  file  for  example  problem  2  is  shown  below: 


1 

1 

1 

0.097 

0 

1 

1 

2 

0.097 

0 

1 

1 

3 

0.097 

0 

1 

1 

4 

0.097 

0 

1 

1 

5 

0.097 

0 

1 

1 

6 

0.097 

0 

1 

1 

7 

0.971 

0 

1 

1 

8 

0.971 

0 

1 

1 

9 

0.971 

0 

1 

1 

10 

0.971 

0 

1 

1 

11 

0.971 

0 

1 

1 

12 

0.971 

0 

1 

1 

13 

-1.029 

0 

1 

1 

14 

-1.029 

0 

1 

1 

15 

-1.029 

0 

1 

1 

16 

-1.029 

0 

1 

1 

17 

-1.029 

0 

1 

1 

18 

-1.029 

0 
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As  in  example  problem  1,  the  command 


vel,  tc2.vel 

and  the  velocity  input  file  can  be  replaced  by  the  lines, 

vel(l:2,l,l)=(. 097,0.) 
vel(3:6,l,l)=(. 971,0.) 
vel (7:8, 1,1)=(. 097,0.) 
vel(9:12,l,l)=(-l. 029,0.) 

The  following  pages  show  the  output  that  would  be  obtained  in  the  .out  file  for  this 
problem. 
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***  FULL  VELOCITY  SYMMETRY  *** 


Example  Problem  3 

The  geometry  for  this  problem  is  the  same  as  for  example  problem  1.  This  example 
illustrates  the  use  of  a  point  source  check.  A  point  source  check  can  be  used  to  determine 
the  adequacy  of  the  surface  discretization  as  well  as  the  need  for  and  the  placement  of 
interior  points.  The  basic  idea  of  a  point  source  check  is  to  consider  the  boundary  surface  as 
an  imaginary  surface  in  an  infinite  fluid  medium  containing  a  finite  number  of  point  sources. 
The  surface  pressures  are  calculated  using  the  CHIEF  formulation  with  a  surface  velocity 
distribution  produced  by  the  point  sources  and  compared  with  the  exact  pressures  due  to  the 
point  sources.  If  the  point  sources  are  inside  the  body,  then  the  CHIEF  computation  is  like 
a  radiation  problem.  If  the  point  sources  are  outside  the  body,  then  the  CHIEF  computation 
is  like  a  scattering  problem.  For  this  example,  we  will  use  one  point  source  located  at  the 
origin,  having  a  complex  weight  of  1+iO.  (This  point  is  outside  the  body.)  For  thin  bodies 
like  this  free-flooded  cylinder,  it  is  usually  better  to  place  the  point  source  outside  the  body, 
since  point  sources  inside  the  body  are  so  close  to  the  surface  that  they  will  produce  a  surface 
velocity  distribution  that  varies  more  rapidly  than  the  ones  likely  to  be  used  in  practice.  The 
input  file  for  this  example  is  shown  below: 


freql=1500/(2*pi) 

nblks=12 

ai=l .68 

ao=2 . 04 

ht=l .2 

htd2= . 5*ht 

rotlim=pi/nblks 

end*params 

rotational ,nblks 

circ_planar,  htd2 
ai , ao , -rotlim , rotlim 

2.1.4.16 
3, 0,0,0 

circ_cylinder,  ao 
-htd2 ,htd2 , -rotlim, rotlim 

4.1.4.16 
3, 0,0,0 

circ_planar,  -htd2 
ai , ao , -rotlim, rotlim 

2.1.4.16 
-3, 0,0,0 

circ_cylinder,  ai 
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-htd2 ,htd2 , -rotlim, rotlim 

4,1,4,16 

-3, 0,0,0 

end*geom 

freq,  freql,  freql,  0. 
point  source  check,  1,1 
0,0, 0,1, 0,1 

The  following  pages  show  the  output  that  would  be  obtained  in  the  .out  file  for  this 
problem. 
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Example  Problem  4 

The  geometry  for  this  example  problem  is  the  same  as  for  example  problem  1.  In  this 
example,  the  radiation  impedance  matrix  relating  surface  forces  and  velocities  and  the  field 
matrix  relating  far-field  pressures  and  surface  velocities  are  computed.  It  is  assumed  that 
the  velocities  have  the  same  symmetry  as  the  geometry.  The  results  are  to  be  written  to  an 
output  file.  The  input  file  for  this  example  problem  is  given  below: 


f reql=1500/ (2*pi) 

nblks=12 

ai=l . 68 

ao=2 . 04 

ht=l . 2 

htd2=.5*ht 

rotlim=pi/nblks 

end*params 

rotational ,nblks 

circ_planar,  htd2 
ai , ao , -rotlim, rotlim 

2.1.4.16 
3, 0,0,0 

circ_cylinder ,  ao 

-ht d2 , ht d2 , -r ot 1 im , r ot 1 im 

4.1.4.16 
3, 0,0,0 

circ_planar,  -htd2 
ai , ao , -rotlim, rotlim 

2.1.4.16 
-3, 0,0,0 

circ_cylinder ,  ai 

-htd2 , htd2 , -rotlim , rotlim 

4.1.4.16 
-3, 0,0,0 

end*geom 

freq,  freql,  freql,  0. 
impedance  matrix, 1,1 
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farfield  pressures, 1,1, 19 
0,90,5,0,0,0 


Example  Problem  5 

In  this  example  problem,  we  consider  a  uniformly  pulsating  sphere  of  radius  1  m  with  a 
radial  velocity  of  1+iO  m/sec.  Figure  8  shows  the  geometry  and  breakup  for  this  example 
problem.  This  example  problem  shows  the  need  and  use  of  interior  points.  We  consider  two 
frequencies  corresponding  to  ka  =  7t/3  and  ka  =  it.  The  value  ka  =  n  corresponds  to  an 
interior  resonance  where  interior  points  are  needed.  The  value  ka  =  7t/3  is  well  below  the 
first  interior  resonance  and,  thus,  does  not  require  interior  points.  Let  us  first  calculate  the 
far-field  pressures  for  these  two  frequencies  with  no  interior  points. 


x3 


Figure  8.  Geometry  for  example  problem  5. 


The  input  file  for  this  calculation  is  shown  below: 

kal=pi/3 
ka2=pi 
a=l .  0 
c=1500 

fu=ka2*c/(2*pi) 
fl=fu/3 
f inc=fu-fl 
uu=pi/6 
end*params 


rot,  6 

sphere ,1.0 
0,pi,-uu,uu 
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12,1,4,16 
3, 0,0,0 
end*geom 


frequency , f 1 , f u , f inc 

!  interior  points, 3 
!  0,0, -.666667 
!  0,0,0 
!  0,0,-666667 

rhs  symmetry 

vel(l:12,l,l)=1.0 

!  or  read  velocities  from  a  file  with:  velocity, tc5.vel 

surface  pressures 

fax-field  pressures, 1, 1,4 
0,180,30,0,0,0 

The  following  pages  show  the  output  that  would  be  obtained  in  the  .out  file  for  this 
problem. 
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0.00  120.00  0. 195511E+06  -.187096E+08  0.187106E+08  -89.40  0.00 
0.00  150.00  0. 195500E+06  -.187095E+08  0.187106E+08  -89.40  0.00 
0.00  180.00  0 . 195495E+06  -.187094E+08  0.187104E+08  -89.40  0.00 


Looking  at  the  directivity  index  (DI)  corresponding  to  ka  =  n  shows  this  answer  is 
incorrect.  Let  us  now  rerun  the  same  problem  with  three  interior  points  located  at  (0,0,- 
2/3),  (0,0,0),  and  (0,0, 2/3).  The  input  file  for  this  calculation  is  the  same  as  the  previous 
one,  with  the  interior  point  control  lines  uncommented.  The  following  pages  show  the  output 
that  would  be  obtained  in  the  .out  file  for  this  example  problem.  It  can  be  easily  verified 
that  the  answers  for  this  case  are  correct. 
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DIRECTION  FOR  PATTERN  NORMALIZATION 


Example  Problem  6 

In  this  example  problem,  we  consider  the  plane-wave  scattering  from  a  rigid  sphere  of 
radius  5  m  and  ka  =  5,  using  two  methods  of  geometry  generation.  Both  methods  use  6-th 
order  rotational  symmetry  with  20  subdivisions  per  symmetry  block.  The  first  method  gen¬ 
erates  the  geometry  using  the  sphere  equations.  The  second  method  generates  the  geometry 
using  axisymmetric  finite  element  inputs.  Figure  9  shows  the  target  strength  versus  angle  for 
both  methods  of  geometry  generation  along  with  the  results  of  an  analytic  series  solution. 


ANALYTIC  CHIEF_SPH  CHIEF_FE 


Figure  9.  Target  strength  versus  angle  for  a  sphere  ( ka  =  5). 

The  contents  of  the  input  file  for  the  spherical  geometry  input  is  shown  below: 


nblks=6 

c=1500 

a=5 

rka=5 

f req=rka*c/ (2*pi*a) 
rotlim=pi/nblks 

end*params 

rot , nblks 
sphere , a 

0 , pi , -r ot 1 im , rot 1 im 

20,1,4,16 

3,0, 0,0 
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end*geom 


frequency , f req , f req , 0 
surface  pressures 
plane  wave, 1,1, 0,0 
far-field,l,3 
0,180,2,0,0,0 


The  input  file  for  the  axisymmetric  finite  element  input  is  exactly  the  same  as  above 
except  that  the  geometry  section  is  replaced  by 

rot , nblks 

fE.axisymmetric,  tc6.crd,  tc6.elt,  4,  16 
3 

end*geom 


The  contents  of  the  nodal  coordinate  file  tc6.crd  are  shown  below. 


1 

0 . 0000E+00 

0.5000E+01 

2 

0 . 3923E+00 

0.4985E+01 

3 

0 . 7822E+00 

0.4938E+01 

4 

0 . 1167E+01 

0.4862E+01 

5 

0 . 1545E+01 

0.4755E+01 

6 

0 . 1913E+01 

0.4619E+01 

7 

0.2270E+01 

0.4455E+01 

8 

0 . 2612E+01 

0.4263E+01 

9 

0 . 2939E+01 

0.4045E+01 

10 

0.3247E+01 

0.3802E+01 

11 

0.3536E+01 

0.3536E+01 

12 

0.3802E+01 

0.3247E+01 

13 

0.4045E+01 

0.2939E+01 

14 

0.4263E+01 

0.2612E+01 

15 

0.4455E+01 

0.2270E+01 

16 

0.4619E+01 

0. 1913E+01 

17 

0.4755E+01 

0. 1545E+01 

18 

0.4862E+01 

0. 1167E+01 

19 

0.4938E+01 

0.7822E+00 

20 

0.4985E+01 

0 . 3923E+00 

21 

0.5000E+01 

-.2186E-06 

22 

0.4985E+01 

- . 3923E+00 

23 

0.4938E+01 

- . 7822E+00 

24 

0.4862E+01 

1167E+01 

25 

0.4755E+01 

1545E+01 
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26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


0.4619E+01 
0.4455E+01 
0.4263E+01 
0.4045E+01 
0.3802E+01 
0.3536E+01 
0.3247E+01 
0.2939E+01 
0.2612E+01 
0.2270E+01 
0 . 1913E+01 
0 . 1545E+01 
0 . 1167E+01 
0 . 7822E+00 
0 . 3923E+00 
-.4371E-06 


■.  1913E+01 
■.2270E+01 
•.2612E+01 
-.2939E+01 
•.3247E+01 
-.3536E+01 
■.3802E+01 
-.4045E+01 
-.4263E+01 
•.4455E+01 
-.4619E+01 
-.4755E+01 
-.4862E+01 
-.4938E+01 
-.4985E+01 
-.5000E+01 


The  contents  of  the  element  file  tc6.elt  is  shown  below. 


1 
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3 
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5 

5 

6 

7 

7 

8 

9 
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18 
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19 

20 

21 

21 

22 

23 

23 

24 

25 

25 

26 

27 
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28 

29 

29 

30 

31 

31 

32 

33 

33 

34 

35 

35 

36 

37 

37 

38 

39 

39 

40 

41 
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Example  Problem  7 

In  this  example  problem,  we  illustrate  the  plane  infinite  rigid  baffle  option.  The  problem 
is  to  compute  the  radiation  impedance  of  a  half-wavelength  diameter  circular  piston  in  a 
plane  infinite  rigid  baffle.  The  impedance  is  to  be  calculated  in  CHIEF  using  four  blocks  of 
rotational  symmetry  with  four  radial  subdivisions  per  block.  The  analytic  solution  for  this 
problem  gave  the  following  result: 

z  =  .952057  x  107  +  *.777030  x  107. 

To  calculate  the  impedance  with  CHIEF  we  use  a  uniform  velocity  distribution  of  1+iO.  The 
impedance  is  then  equal  to  the  complex  power  given  in  the  .out  file.  The  value  computed 
with  this  breakup  is 

z'=  .957883  x  107  +  *.780259  x  107, 

which  agrees  with  the  analytic  solution  to  within  about  0.5%.  The  input  file  for  this  run  is 
shown  below: 


freq=1500/(2*pi) 

nblks=4 

a=pi/2 

rotlim=pi/nblks 

end*params 

rot ,nblks 
circ_planar,0 
0 , a , -rotlim,rotlim 
4, 1,2,6 
3 

end*geom 

frequency , f req, f req, 0 

inf  baffle 

rhs  symmetry 

velocity , t c7 . vel 
!  or  vel(l:4,l,l)=1.0 

surf ace, 3 
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The  velocity  file  tc7.vel  is  as  follows 


1  1  1  1.  0. 

1  1  2  1.  0. 

1  1  3  1.  0. 

1  1  4  1.  0. 
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Example  Problem  8 

This  example  problem  shows  how  to  compute  near-field  pressures.  The  radiating  surface 
consists  of  a  circular  cylinder  with  flat  circular  end  pieces.  The  cylinder  radius  is  1  m  and  the 
height  is  2  m.  We  choose  the  origin  at  the  center  of  the  cylinder  and  let  the  z-axis  coincide 
with  the  cylinder  axis.  The  end  surface  at  2  =  1  is  given  a  velocity  of  1+iO  and  the  rest  of 
the  surface  has  zero  velocity.  We  want  to  compute  the  pressure  at  the  two  points  (0,0,1-25) 
and  (0,0, 1.5)  on  the  axis  near  the  piston.  The  input  file  for  this  problem  is  shown  below: 

freql=100 

frequ=1500 

freqinc=1400 

nblks=12 

ao=l .0 

ht=2 . 0 

htd2= . 5*ht 

rotlim=pi/nblks 

end*params 

rotational ,nblks 

circ_planar,  htd2 
0 , ao , -rot 1 im , rot 1 im 

16.1.4.16 
3 

circ_cylinder,  ao 

-htd2 , htd2 , -r ot lim , rot 1 im 

32.1.4.16 
3 

circ_planar,  -htd2 
0,ao,-rotlim,rotlim 

16.1.4.16 
-3,0, 0,0 

end*geom 

freq,  freql,  frequ,  freqinc 
rhs  symmetry 
vel(l:16,l,l)=1.0 
surface  pressures 
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neaxfield  pressures, 2, 3 

0,0,1.25 

0,0, 1.5 

Note  that  you  only  have  to  input  the  nonzero  velocities.  The  results  from  the  file  tc8.out 
are  shown  on  the  following  page 
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. 163006E+07  -. 110086E+07  0.196698E+07  -145.97  0.00 

. 267511E+07  0 . 117691E+07  0.292256E+07  156.25  3.44 


DIRECTION  FOR  PATTERN  NORMALIZATION 
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Example  Problem  9 

This  example  problem  is  like  example  problem  2  except  that  both  surface  and  far-field  A 
and  B  matrices  are  output.  In  this  case,  it  is  not  necessary  to  input  a  velocity  distribution. 
The  control  section  of  the  input  file  for  this  example  is  shown  below: 


freq,  freql,  freql,  0. 

a_b  matrices 

surface  matrices 

farfield  matrices, 1,0 
0,90,5,0,0,0 
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APPENDIX:  BASIC  CHIEF  THEORY 


The  principal  quantity  of  interest  in  classical  (linear)  acoustics  is  the  acoustic  pressure 
p(x,  t )  which  is  a  function  of  both  position  and  time.  The  acoustic  pressure  is  the  excess 
pressure  over  the  hydrostatic  pressure  due  to  compressional  waves  in  the  medium.  The 
acoustic  pressure  p(x,  t)  satisfies  the  wave  equation, 


A  p  — 


1  d2p 
c 2  dt2 


(1) 


Here,  A  is  the  spatial  Laplacian  operator  and  c  is  the  speed  of  sound  in  the  fluid  medium.  In 
acoustics,  it  is  often  more  convenient  to  work  in  the  frequency  domain  than  in  the  time  do¬ 
main.  By  frequency  domain  we  mean  the  Fourier  transform  domain.  The  Fourier  transform 
P  of  p  is  defined  by  the  pair  of  relations, 

/OO 

P(x,u))etult  df  (inverse  Fourier  transform) ,  (2a) 

-OO 

where 


P(x,u>)  =  /  p(x,  t)e  iwt  dt  (Fourier  transform).  (2b) 

J  —  OO 

The  transform  variable  u)  is  called  the  angular  frequency,  and  the  related  quantity  /  =  u)/ (27 r) 
is  called  the  frequency.  Note  that  whereas  the  pressure  p  is  a  real  quantity,  the  transformed 
pressure  P  is,  in  general,  complex.  The  Fourier  transformed  pressure  P(x,  u)  satisfies  the 
reduced  wave  equation, 

A  P  +  k2P  =  0,  (3) 

where  k  =  uj/c  is  the  acoustic  wavenumber. 

In  scattering  problems,  the  incident  pressure  pinc  is  the  pressure  that  would  be  present 
if  the  scattering  object  were  removed.  If  the  source  of  the  incident  pressure  is  far  from  the 
scattering  object,  pmc  is  usually  taken  to  be  a  plane  wave  since  the  wave  is  approximately 
planar  in  the  vicinity  of  the  scatterer.  If  p  is  the  total  acoustic  pressure  present  in  a  scattering 
problem,  then  the  scattered  pressure  ps  is  defined  by 


p=p-p 


(4) 


Clearly,  the  incident  and  scattered  pressures  also  satisfy  the  acoustic  wave  equation. 

To  uniquely  determine  the  pressure,  it  is  necessary  to  specify  some  initial  and  boundary 
conditions  in  addition  to  the  requirement  that  the  wave  equation  be  satisfied.  In  the  time 
domain,  we  usually  specify  initial  conditions  throughout  the  region  of  interest  for  both  p 
and  p  as  well  as  boundary  conditions  on  either  p  or  the  normal  velocity  v.  In  the  time 
domain,  these  boundary  and  initial  conditions  are  sufficient  to  uniquely  determine  p  for 
exterior  (unbounded  region)  problems.  In  the  frequency  domain,  we  usually  specify  either 
the  pressure  transform  P  or  the  normal  velocity  transform  V  on  the  boundary.  For  exterior 
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problems  in  the  frequency  domain,  it  is  necessary  to  supplement  the  boundary  conditions 
with  a  radiation  condition  at  infinity  to  eliminate  the  possibility  of  waves  coming  in  from 
infinity.  This  condition  usually  takes  the  form, 

}™r(^+ikP)=0’  (5) 

due  to  Sommerfeld.  In  effect,  this  condition  says  that  the  pressure  field  looks  locally  like 
an  outgoing  spherical  wave  at  large  distances  from  the  sound  source.  This  condition  can  be 
eliminated  if  the  fluid  is  assumed  slightly  lossy  (complex  sound  velocity).  If  the  incident  field 
in  a  scattering  problem  is  considered  a  plane  wave,  then  neither  the  incident  pressure  nor  the 
total  pressure  satisfies  the  outgoing  radiation  condition.  Therefore,  in  scattering  problems, 
the  outgoing  radiation  condition  is  applied  to  the  scattered  pressure. 

The  CHIEF  method  for  numerically  solving  acoustic  problems  is  based  on  the  reduction 
of  a  boundary- value  problem  for  the  wave  equation  in  the  frequency  domain  to  an  equivalent 
integral  equation  over  the  boundary  surface.  This  integral  equation  is  obtained  from  the 
Helmholtz  integral  relations: 


Js[-iupG(x,t)V(0  -  P({)^(: t,{)]  dS(() 


— P(x )  x  outside  S 
—\P{x)  x  on  S 
0  x  inside  S 


where  G(x,£)  is  the  free-space  Green’s  function  defined  by 


(6) 


G(*,0 


e— ifc|x— 1| 

47r|z-f|’ 


(7) 


n  is  the  outward  unit  normal  to  S,  and  \x— £|  is  the  distance  between  x  and  These  relations 
come  from  the  divergence  theorem,  the  fact  that  both  P  and  G  satisfy  the  wave  equation,  and 
the  nature  of  the  singularity  of  G.  A  good  discussion  of  these  integral  relations  is  contained 
in  [Baker  and  Copson,  1950].  If  the  normal  velocity  V  is  specified  on  the  boundary  S ,  then 
the  middle  relation  in  equation  (6)  provides  an  integral  equation  for  the  surface  pressure  P, 
for  example, 


HO  -  JsP(0^(C0dS(0  =  iojpJsG(COV(OdS(0  C  on  5.  (8) 

Once  this  integral  equation  is  solved  for  the  surface  pressure,  the  pressure  at  any  point  in 
the  exterior  region  can  be  obtained  using  the  first  relation  in  equation  (6).  However,  for 
the  exterior  problem,  there  is  an  infinite  discrete  set  of  frequencies  /1;  /2, . . .  for  which  the 
integral  equation  does  not  have  a  unique  solution.  It  turns  out  that  these  frequencies  are  the 
resonance  frequencies  of  the  interior  problem  for  the  boundary  condition  P  =  0  on  S.  As  one 
might  expect,  numerical  methods  based  on  these  integral  equations  have  a  problem  when  the 
frequency  is  close  to  one  of  these  interior  resonance  frequencies.  There  is  a  large  volume  of 
literature  describing  various  methods  for  overcoming  this  difficulty  [Benthien  and  Schenck, 
1997].  The  resonance  frequencies  /i,  f2, . . .  are  not  resonances  of  the  problem  of  interest. 
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There  are  no  resonances  of  the  exterior  problem.  In  fact,  the  exterior  problem  has  a  unique 
solution  at  all  frequencies.  The  interior  resonance  problem  is  a  mathematical  difficulty,  not 
a  physical  one.  What  happens  at  these  frequencies  is  that  the  integral  equation  is  no  longer 
equivalent  to  the  original  exterior  acoustic  boundary  value  problem. 


CHIEF  uses  a  piecewise  constant  approximation  to  the  integral  equation  (8).  In  this 
approximation,  the  boundary  surface  S  is  subdivided  into  a  finite  number  of  subareas 
Si,  S2, ... ,  Sn •  On  each  subarea  Sn  the  pressure  is  approximated  by  a  single  value  Pn  and  the 
normal  velocity  is  approximated  by  a  single  value  Vn.  With  this  approximation  the  integral 
equation  (8)  becomes 


ip(C)-X>„/  7p(C.f)<tf(0  =  »"p2>"  /  G(U)dS((). 

n—1  '*Sn  U  n=l  ^Sn 


(9) 


Evaluating  equation  (9)  at  a  set  of  reference  points,  Ci>  C2 , . . . ,  Cat  (one  on  each  subarea),  we 
obtain  the  system  of  equations 


N 


N 


^  ^  AmnPn  ^  ^  BmnVn 


m 


1,2,. ..,7V, 


71=  1 


71=1 


(10) 


where 

A  mn  =  \5mn-  Js  ^(UO  dS{(), 
Bmn  =  iu,p[  G(U,()dS((), 

JSn 

and  5mn  is  the  Kronecker  delta  defined  by 


(11) 

(12) 


_ _ m*n  (13) 

II  m  =  n. 

is  taken  at  the  center  of  the  m-th  subarea.  The  system  of  equations  (10)  can  be  written 
in  the  matrix  form, 

AP  =  BV,  (14) 

where  P  is  an  N-vector  whose  n-th  component  is  Pn  and  V  is  an  N-vector  whose  n-th 
component  is  Vn.  If  the  velocities  Vi,  V2, . . . ,  Vn  are  prescribed,  then  the  system  of  equations 
(10)  can  be  solved  for  the  surface  pressures  Pi,P2,...,  Pn- 


The  expression  in  equation  (6)  for  the  pressure  at  an  arbitrary  field  point  can  be  approx¬ 
imated  using  the  same  discrete  approximations,  yielding 


N 


P(x)  =  Y,[A„(x)Pn  +  BJx)V,t\, 


(15) 


71=1 


where 


A„(x)  =  J  and  Bn(x)  =  iup G(x,()dS(t).  (16) 
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This  expression  can  be  used  to  calculate  the  pressure  at  an  arbitrary  field  point  x  once  the 
surface  pressures  and  normal  velocities  have  been  determined.  Evaluating  equation  (15)  at 
a  set  of  field  points  aq, . . . ,  we  obtain 

N 

P(%m)  =  ^  ^[An(Xm)Pn  d*  Bn(xm)V^]  7Tb  =  1,  .  .  .  ,  M.  (17) 

n=l 

Equation  (17)  can  be  written  in  the  matrix  form, 

pfld  =  Amp  +  pfld  y-  (i8) 

where  Pfld  is  the  vector  of  field  pressures,  P  is  the  vector  of  surface  pressures,  V  is  the  vector 
of  normal  surface  velocities,  and  the  field  matrices  AM,  BM  are  given  by 

=  An(xm)  (19) 

B"  =  Bn(xm).  (20) 

CHIEF  also  contains  far-field  approximations  for  these  field  equations  when  the  field  points 
are  far  from  the  radiating  surface.  Equation  (14)  can  be  solved  for  P  to  obtain  the  impedance 
relation, 

F  =  DP  =  ZV,  (21) 

where  D  is  the  diagonal  matrix  of  surface  subdivision  areas  and  Z  is  the  impedance  matrix 
given  by 

Z  =  DA~XB.  (22) 

Substituting  equation  (21)  into  equation  (18),  we  obtain 

pfld  =  +  pfldjy;  (23) 

The  matrix  on  the  right-hand  side  of  equation  (23)  is  sometimes  referred  to  as  the  field 
matrix. 

In  the  vicinity  of  one  of  the  frequencies  for  which  the  underlying  integral  equation  does 
not  have  a  unique  solution,  the  system  of  equations  (10)  becomes  numerically  ill-conditioned. 
One  method  of  overcoming  this  difficulty  is  to  add  some  additional  equations  obtained  by 
evaluating  a  discretized  version  of  the  last  relation  in  equation  (6)  at  a  selected  number  of 
interior  points.  The  overdetermined  system  of  equations  is  then  solved  in  a  least  squares 
sense.  A  more  complete  discussion  of  this  approach  to  solving  acoustic  problems  is  contained 
in  [Benthien,  Barach,  and  Gillette,  1988]. 

While  it  is  important  to  maintain  the  capability  of  handling  quite  general  surfaces,  it  is 
also  important  to  realize  that  significant  reductions  in  computation  time  can  be  achieved  by 
taking  advantage  of  any  symmetries  that  might  be  present  in  the  problem.  In  the  CHIEF 
approach,  symmetry  manifests  itself  in  the  structure  of  the  A  and  B  matrices.  The  previously 
referenced  CHIEF  manual  shows  that  each  type  of  symmetry  results  in  the  A  and  B  matrices 
commuting  with  the  appropriate  symmetry  operator  E,  for  example, 

AH  =  E  A  and  BE  =  E  B.  (24) 
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It  follows  from  these  commutation  relations  that  the  eigenspaces  of  £  are  invariant  under 
both  A  and  B.  Thus,  A  and  B  become  block  diagonal  when  expressed  in  terms  of  a  basis 
consisting  of  eigenvectors  of  £.  The  sizes  of  the  diagonal  blocks  are  the  dimensions  of  the 
eigenspaces  of  E. 

As  an  example  of  how  symmetry  is  treated  in  CHIEF,  consider  the  special  case  of  one 
plane  of  symmetry.  The  symmetry  operator  for  one  plane  of  symmetry  is 


E  = 


(25) 


assuming  that  symmetric  elements  are  numbered  in  the  same  order.  The  fact  that  A  com¬ 
mutes  with  E  implies  that  A  has  the  block  form, 


A  = 


A\  A2 

i  A-2  Ai 


(26) 


The  matrix  X  of  eigenvectors  of  E  is  given  by 


I  I 
I  -I 


(27) 


The  inverse  X  1  of  X  turns  out  to  be  ^ X .  It  is  easily  verified  that  X  '  AX  and  X  ^BX 
have  the  block  diagonal  form, 


X^AX  = 


Ai  +  A2  0 
0  Ai  —  A2 


,  X~lBX  = 


B\  +  B2 
0 


0 

B\  —  B2 


Thus,  the  equation  AP  =  BV  can  be  written  as 

(X~1AX)P  =  (X-'BXW, 

where 

P  =  X-'P  and  V  =  X~ 1V. 
Equation  (29)  can  be  written  in  the  partitioned  form, 


Ai  +  A2  0  \  (PA  I  B\  +  B2 

0  Ai  —  A2  J  \p2  /  \  0 


0 


b1-b2  \v2 


(28) 

(29) 

(30) 

(31) 


Thus,  the  original  system  of  equations  reduces  to  the  pair  of  smaller  systems  of  equations: 

(A1+A2)Pi  =  (B1  +  B2)V1  (32) 

C A: l  -  A2)P2  =  ( B x  -  B2)V2.  (33) 
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Once  these  equations  are  solved,  the  solution  of  the  original  system  of  equations  is  obtained 
by 

Pi  =  Pi  +  P2  (34) 

P2  =  A  -  P2.  (35) 

Since  solution  time  goes  roughly  as  the  cube  of  the  number  of  equations,  it  is  faster  by 
a  factor  of  four  to  solve  the  two  smaller  systems  rather  than  the  one  larger  system.  The 
symmetry  can  also  be  used  to  reduce  the  time  needed  to  generate  the  matrices.  The  CHIEF 
program  can  take  advantage  of  one,  two,  or  three  planes  of  symmetry  as  well  as  any  finite 
order  of  rotational  symmetry  in  the  surface.  CHIEF  can  also  take  advantage  of  symmetry 
in  the  velocity  distribution  when  it  exists.  Table  A-l  shows  the  time-reduction  factors  for 
various  types  of  symmetry. 


Table  A-l.  Time  reduction  factors  in  CHIEF  for  various  types  of  symmetry. 


Symmetry  Type 

Solution 

Using  velocity  symmetry 

Matrix  generation 

One  plane  of  symmetry 

4 

8 

2 

Two  planes  of  symmetry 

16 

64 

4 

Three  planes  of  symmetry 

64 

512 

8 

A-fold  rotational  symmetry 

N2 

N3 

N 

The  reader  should  see  the  appendix  to  the  original  CHIEF  manual  for  more  details  on 
the  theory  behind  the  CHIEF  program. 
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